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PREFACE 
The  application  of  artificial  earth  satellites  in  geodesy  can  be 
either  dynamic  or  geometric o 

In  the  dynamic  approach^  an  evaluation  of  the  perturbations  in 
the  orbital  motion  of  a  satellite  can  furnish  data  from  which  the  coef- 
ficients of  the  spherical  harmonic  expansion  of  the  potential  function 
can  be  determinedo  Using  this  approach,,  valuable  information  on  the 
gravitational  field  of  the  earth  can  be  obtained » 

In  the  geometric  treatment^  the  artificial  earth  satellite  is 
used  as  a  reference  point  or  a  triangulation  station  in  the  sky0  In 
this  way  the  positions  of  observation  stations  on  the  surface  of  the 
earth  can  be  determinedo  Some  aspects  of  the  geometric  approach  will 
be  discussed  in  this  paper© 

The  purpose  of  this  thesis  is  to  serve  as  a  demonstration  of  the 
material  studied  during  its  preparation »  I  would  like  to  emphasize  that 
the  reader  should  not  construe  the  methods  and  procedures  used  herein 
as  being  of  a  standard  type  because  of  the  following  reasons? 
lo  The  computed  orbit  from  three  observations  at  a 
given  station  at  short  intervals  is  assumed  to  be  correct 
and  accurate 5 

2o  It  is  assumed  that  the  observed  coordinates  of  the 
satellite  taken  from  Smithsonian  Institution  Astrophysical 
Observatory  Special  Report  No,  66  and  the  computed  orbital 
elements  are  referred  to  the  same  mean  coordinate  system^ 

3o  The  Gaussian  Method  used  to  determine  the  orbit  is  an 
approximate  method| 

k»     It  is  assumed  that  only  two  approximations  are  necessary 
in  the  differential  process  of  computing  the  geographic  latitude 
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and  longitude • 
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1.   INTRODUCTION 

Since  the  early  days  of  history  man  has  shown  a  keen  interest  and 
an  investigative  concern  about  the  earthy  its  characteristics!,  and  its 
dimensions o  Since  the  third  century  B»C0  when  Pythagoras  stated  that 
the  earth  was  of  spherical  shape ,  many  attempts  have  been  made  to 
measure  its  size  and  shape  and  to  reduce  the  physical  parameters  and 
features  on  to  maps  and  chartso  Of  the  many  existing  earth  sciences., 
geodesy  is  the  one  that  primarily  contains  this  field* 

The  inception  of  geodetic  science  can  be  traced  back  to  Eratosthenes 
(circa  2?6  -  2Sk  BoCo)  who  was  the  first  man  to  measure  the  size  of  the 
eartho  His  method  of  measuring  the  length  of  a  meridianal  arc  by  taming 
the  travel  of  a  camel  caravan  between  Alexandria  and  Syene  (Aswan)  and 
measuring  the  angle  of  incidence  of  the  rays  of  the  sun  shining  into  a 
well  at  these  places  was  crude  but  effective  „  The  degree  of  accuracy 
he  obtained  was  remarkable  when  we  consider  the  tools  he  had  available  <, 

The  science  showed  little  development  between  this  time  and  the 
seventeenth  century  <>  At  this  time  the  advent  of  the  telescope,  loga~ 
rithmic  tables  and  the  method  of  triangulation  (developed  by  the  Dutch 
scientist  Snellius  in  l6lS>)  took  place  and  geodesy  was  reborn » 

The  history  of  geodesy  can  be  divided  into  periods  based  on  the 
figure  of  the  earth  that  was  used  as  a  reference  during  the  respective 
periods o  The  studies  and  experiments  of  Pythagoras 9   Aristotle s   and 
Eratosthenes  brought  the  "flat%  or  plane  era  to  an  end  and  introduced 
the  spherical  era«  Newton  and  Huygens9  in  the  eighteenth  century^ 
introduced  the  ellipsoidal  period  which  fin  ally  phased  into  the 
geoidal  and  telluroidal  periods  of  this  century  o 
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With  the  launching. of  the  first  artifici&L  satellite  and  the  begin- 
ning of  the  space  age  in  this  decade  a  new  era  in  the  field  of  geodesy 
is  unfolding o     By  using,  the  geoid  as  a  reference  surface  we  can  determine 
the  shape  of  the  earth,  but  not  its  size0     The  orbiting  artificial  earth 
satellites  will  not  replace  the  geoid  or  the  telluroid  but  will  provide 
the  geodesist  with  a  convenient  reference  station  or  "base-line  in  the 
sky"  that  can  be  used  to  measure  both  the  size  and  shape  of  the  earth© 

Since  the  launching  of  our  first  artificial  satellite  (Explorer 
I)  on  31  January  1958,  the  field  of  space  science  has  expanded  tremen- 
dously and  with  the  successful  results  of  the  early  and  present  projects 
the  limits  in  this  science  are  beyond  imagination*     At  the  present, 
the  projects  already  underway  and  those  planned  for  the  immediate 
future  suggest  the  following  three  areas  of  space  activity?  research, 
application^  and  exploration*      [1,  P<>5> 

The  applications  include  meteorology,,  communications-,  and  geodesy© 
No  new  principles,  yet  to  be  discovered  are  necessary;  no  major  tech- 
nological problems  stand  in  our  way* 

In  the  field  of  meteorology,  the  satellites  can  provide,  through 
telemetry*  information  on  pressure,  temperature,  cloud  coverage  and 
stoaaa  patterns  and  also  assist  in  predicting  such  mankillers  as 
hurricanes,  cyclones,  etc»    The  day  to  day  practical  use  of  weather 
satellites  is  obvious  o 

In  communications^  satellite  contributions  are  numerous: 
telephony,  long-range  radio  communications,  and  international  tele- 
vision o 

This  brings  us  to  the  third  application  of  space  activity: 
geodesy  <>     Geodetic  satellites  can  serve  both  research  and  application 


3 
interests  o  As  to  the  latter,  navigation  and  surveying,  come  closest  to 
the  immediate  interests  of  mankind.  The  geodesist  can  use  the  satellite 
in  both  a  dynamic  and  a  geometric  approach.  Dynamically,  the  pertur~ 
bations  in  the  orbits  of  the  satellites  caused  by  changes  in  the  earth's 
gravitational  attraction  can  be  mathematically  computed  and  analyzed 
and  thus  provide  valuable  information  to  assist  in  determining  the 
gravity  field  of  the  earth0  By  observing  the  changes  in  the  orbits 
of  many  satellites  over  a  long  period  of  time,  information  on  regional 
gravity  anomalies  can  be  obtained ■  From  these  data  we  can  compute 
the  shape  of  our  equipotential  reference  surface,  the  geoido  Further, 
observations  of  the  nodal  motion  or  the  rate  of  change  of  the  plane 
of  the  orbit  of  the  satellite  gives  us  information  on  the  flattening 
of  the  earth.  Therefore,  in  the  dynamic  approach  we  can  determine  the 
shape  of  the  eartho 

The  size  of  the  earth  can  be  realized  in  the  geometric  approach. 
In  this  process  the  satellite  will  serve  as  a  reference  point  or  a 
triangulation  station  in  the  sky„  Thus  we  can  connect  triangulation 
nets  or  stations  that  are  widely  separated  by  oceans  and  continents 
or  are  located  in  inaccessible  areas  where  the  conventional  methods  of 
triangulation  would  be  awkward  and  cumbersome «  If  so  desired,  the 
existing  geodetic  systems  could  be  reduced  to  a  common  world  geodetic 
network  or  datum. 

In  this  thesis  a  phase  of  the  geometric  use  of  the  artificial 
earth  satellite  in  geodesy  will  be  discussed.  In  the  geometric 
approach  the  geodesist  can  use  the  satellite  in  two  basic  ways.  One 
method  involves  the  observation  of  the  satellite  simultaneously  from 
two  or  more  stations.  The  other  method  requires  that  the  orbit  of  the 
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satellite  be  precisely  known <>     Then  knowing  the  orbital  parameters  or 
characteristics  and  predicting  their  values  for  a  desired  time  we 
observe  the  satellite  from  any  number  of  stations  on  the  eartho 

The  discussion  of  the  problem  "will  be  treated  in  the  following 
ways 

First5  the  orbit  of  the  satellite  Explorer  I  (19£8  Alpha)  *ill  be 
computed  from  the  right  ascension  and  declination  of  the  satellite  as 
observed  at  three  different  times  from  the  tracking  station  at  Olifants- 
fontein^  South  Africa..     The  three  observations  have  been  made  at  very 
short  intervals  during  the  same  revolution  of  the  satellite  <>     The 
values  of  right  ascension  and  declination  are  nppttmrmt  ones  « 
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from  camera  films  of  Baker=»Nunn  tracking  cameras o     The  geographic 
coordinates^  latitude  (  m )  and  longitude  {  \  )  of  the  tracking  station 
are  known  and  its  height  (H)  above  sea  level  is  given,. 

Second^  the  same  satellite  will  be  observed  again  from  a  station 
on  the  island  of  Curacao  in  the  Netherland  West  Indies  in  the  Caribbean0 
The  geographic  coordinates  of  this  station  will  be  assumed  to  be  unknownc 
Then  using  the  computed  orbital  parameters  and  observing  the  satellite's 
right  ascension  and  declination  from  this  station^,  the  geographic  co- 
ordinates of  this  unknown  station  will  be  found o 

Thirds  the  forward  and  back  azimuths  and  the  geodesic  distance 
between  the  two  stations  will  be  computedo 

This  connection  between  the  two  points  is  called  a  tie  and  in 
this  instance s  an  inter-continental  tie<>     The  distance  between  the 
stations  is  approximately  6(!000  nautical  miles o    Tho  procedure  dioouoood 
in  tho  following  chapters  can  be  applied  to  any  satellite  observed  ftw 
any  looation  on  tho  surface  of  the  eartho 


20       DEFINITIONS  AND  ELEMENTS  OF  THE  ORBIT 

The  motion  of  an  artificial  satellite  in  its  orbit  about  the 
earth  follows  the  same  lairs  that  apply  to  the  movement  of  the  planets 
in  relation  to  the  sun  in  the  solar  system© 

Johannes  Kepler  (l£71  -  1630),  a  German  astronomer  and  mathema- 
tician, derived  the  laws  of  planetary  motion  from  the  observations  of 
Tycho  Brahe   (l£li6  -  1601),  Danish  astronomer,  with  whom  he  worked  as 
an  assistant*, 

Kepler's  Laws  ares 

I  The  orbit  of  each  planet  is  an  ellipse  with 
the  sun  at  one  of  its  foci. 

II  The  line   (radius  vector)  joining  the  sun  to 
each  planet  sweeps  over  equal  areas  of  its 
ellipse  in  equal  times • 

III  The  squares  of  the  orbital  periods  of  the 
planets  are  proportional  to  the  cubes  of  their 
mean  distances  from  the  sun  (foci)© 

The  relationship  of  an  artificial  satellite  with  respect  to  the 

earth  is  the  same  as  that  existing  between  the  planets  and  the  sun  as 

> 
stated  in  the  above  laws0 

In  other  words,  we  can  replace  the  term  planet  with  satellite 
and  replace  the  sun  with  the  earth  and  the  laws  are  still  applicable • 
Now,  since  we  know  the  basic  form  of  the  orbit  is  an  ellipse  we  can 
easily  describe  the  characteristics  of  the  orbit .     These  orbital 
characteristics  are  known  as  parameters  and  can  be  considered  to  be 
constant  quantities  for  a  particular  satellite  <,     The  first  two  para- 
meters are  properties  of  the  ellipse  itself  and  these  are  the  semi- 
major  axis,  a,  and  the  eccentricity,  e0     These  two  terms  define  any 
ellipse • 


Figure  1. 


Referring  to  Figure  1.,  the  eccentricity  is  the  ratio  existing 
between  the  length  of  the  semi-major  axis  and  the  distance  of  the 
foci,  c,  from  the  center,  o,  ors 

e  »  c/a 
The  relationship  e^  -  (a^  -  t >  )/a^,  where  b  is  the  minor  semi-axis, 
also  exists.     The  eccentricity  of  a  circle  is  zero  and  that  of  a 
parabola  is  equal  to  one.     The  eccentricity  of  an  ellipse  is  less 
than  one  but  greater  than  zero.     The  eccentricity  of  the  orbit  of 
a  satellite  is  nearly  circular,  that  is,  slightly  greater  than  zero. 

Two  other  parameters  of  the  orbit  are  the  angle  of  inclination, 
i,  and  the  right  ascension  of  the  ascending  node,   St  •       These  two 
values  define  the  plane  of  the  orbit  in  space.     Referring  to  Figure  2, 
the  angle  i  is  measured  from  the  plane  of  the  equator  to  the  plane  of 
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FIGURE  2 
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the  orbit.     It  can  have  a  value  between  0°  and  180°.     The  angle  it    is 
measured  eastward  along  the  equatorial  plane  from  the  vernal  equinox, 
to  the  ascending  node.     It  can  have  any  value  between  0°  and  3^0° • 
The  ascending_  node  is  that  point  where  the  satellite  crosses  the  equa- 
torial plane  from  south  to  north. 

Another  parameter  of  the  orbit  is  the  argument  of  the  perigee, 
w,  or  the  angle  measured  from  the  equatorial  plane  at  the  ascending 
node  to  the  major  semi-axis  (the  extension  of  which  contains  the 
perigee).     The  angle  w  is  measured  in  the  direction  of  the  satellite's 
motion  and  gives  us  the  orientation  of  the  orbit  in  its  plane. 

The  sixth  and  last  parameter  is  the  time,  To.     Either  the  time 
of  passage  of  the  satellite  through  the  ascending  node  or  the  time 
of  passage  through  perigee  may  be  chosen.     The  latter  time  is  usually 
the  one  used  in  most  computations. 

In  summarizing;     the  a  and  e  specify  the  size  and  shape  of  the 
orbit;  i  and    ^  1   specify  the  orientation  of  the  orbital  plane  in  space | 
w  orients  the  orbit  in  its  plane j  T0 furnishes  us  with  the  value  to 
locate  the  satellite  in  its  orbital  plane  relative  to  the  perigee. 

To  position  the  satellite  at  any  specific  time  we  must  resort  to 
a  system  of  angular  values  called  anomalies,  of  which  there  are  threes 
mean  M,  eccentric  E,  and  true  V« 
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Figure  3< 


Refer  to  Figure  3»     Angle  QFN  is  the  true  anomaly  V,  which  is  the 
angular  distance  of  the  satellite  measured  from  the  line  of  apsides 
along  the  plane  of  the  orbit.     Since  the  orbit  of  the  satellite  is 
an  ellipse,  the  radius  vector  FQ,  or  the  distance  of  the  satellite 
from  the  foci,  will  not  be  constant <>     Therefore,  the  angular  velocity 
is  not  constanto     This  must  be  so  if  Kepler's  second  law  is  to  be 
fulfilled*     If  the  angular  velocity  of  the  satellite  was  constant, 
the  radius  vector  would  sweep  at  an  average  rate  and  the  eccentric 
anomaly  would  equal  the  mean  anomaly.     The  mean  anomaly  is  the  theo- 
retical angular  velocity  of  a  satellite  moving  uniformly  through  3^0° 
of  arc  per  orbital  periodo 
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The  symbols  and  definitions  of  terms  used  in  this  paper  are  listed 

be lows 

a  Semi-major  axis  of  the  orbit  of  the  satellite 

(Chapters  h9  £); 
Semi-major  axis  of  the  International  Ellipsoid  (earth) 
(Chapters  1,  2S  3) 

b  Semi-minor  axis  of  the  orbit  of  the  satellite 

(Chapters  h$  5)} 
Semi-minor  axis  of  the  ellipsoid  (Chapters  1*  2,  3) 

c  Distance  from  center  of  orbit  to  foci 

e  Eccentricity  of  orbital  ellipse   (Chapter  li,5)l 

Eccentricity  of  International  Ellipsoid 
(Chapters  la  2,  3) 

f  Flattening  of  the  ellipsoid 

fsg5h  Direction  cosines  (used  together) 

i  Angle  of  inclination  of  the  orbital  plane 

k  Gaussian  gravitational  constant 

n  Mean  daily  motion  of  satellite  (velocity  of  the 

mean  anomaly)  in  revolutions* 

p  Auxiliary  parameter  of  orbital  ellipse  (replaces  a) 

q  Geocentric  distance  of  satellite  at  perigee 

(Chapter  5)  J 
Time  intervals   (Ti  +   T3)k         (Chapter  k) 

r  Distance  of  satellite  from  observer 

s  Geodesic  distance 

t  Time  intervals   (T3  -  T^)k 

u  Auxiliary  parameter  of  orbit 

x,y5z  *     Geocentric  space  rectangular  coordinates 

xa5ya,za  Approximate  space  rectangular  coordinates  of 
the  observer 

xosyOj,zo  Space  rectangular  coordinates  of  the  observer.. 
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A,B,C,D,E,FSG  Various  terms  used  in  Sodano's  Fourth  Method 

E  Eccentric  anomaly  of  satellite 

G  Newton  °s  Universal  gravitational  constant 

H  Elevation  above  sea  level 

K5  Kl  Gaussian  determinants 

Kj  (j  -  1S253*U<)5>56)  Constants  for  International 

Ellipsoid  used  in  SodanCs  Fourth  Method,, 

L  Difference  in  longitude  between  two  points 

on  the  ellipsoid 

M  Mean  anomaly  of  satellite 

N  Radius  of  curvature  in  the  prime  vertical 

R  Geocentric  distance  of  satellite 

Tj  Time  of  observation       (5  =  1,  2^3) 

T0  Time  of  Epoch 

V  True  anomaly  of  satellite 

X,YaZ  Space  rectangular  coordinates  of  the  satellite « 

(X  Flattening  of  the  ellipsoid  (Chapter  6) 

(X  Geocentric  Right  Ascer&on  -  the  angular  distance 

at  the  center  of  the  earth  measured  eastward 
along  the  equator  from  the  apparent  vernal 
equinox  to  the  celestial  meridian  passing 
through  the  satellite • 

i 
(X  Apparent  Right  Ascension  -  measured  at  the 

position  of  the  observer. 

ry  Apparent  right  ascension  obtained  through 

c  the  approximate  coordinates  of  observer© 

(Xio  Forward  Azimuth 

G(p/  Back  Azimuth 

R  Parametric  latitude  (reduced  latitude) 

R  Maximum  latitude  of  the  geodesic  o 

Y  Normal  gravity., 
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£  Geocentric  declination  -  the  angular  distance  at 

^  the  center  of  the  earth  measured  northward  or 

southward  from  the  equator  to  the  satellite.     A 
satellite  below  the  equator  would  have  a  south 
declination  and  the  angle  would  have  a  minus  value o 

£  Apparent  declination  -  measured  at  the  observer's 

^  position. 

£  Computed  apparent  declination  obtained  through 

^c  the  approximate  coordinates  of  the  observer. 

Q  Spherical  distance 

Q  True  spherical  distance 

Qa,  Greenwich  hour  angle  of  the  apparent  vernal  equinox, 

C/fl  Local  hour  angle  of  the  apparent  vernal  equinox 0 

1  Geographic  longitude  (positive  westward  and 

negative  eastward) 

P  Geocentric  distance  of  observer 

rf  The  angular  motion  of  the  satellite  during  the 

interval  between  observations.  Usually  expressed 
in  radians. 

(f\  Geographic  latitude  (north  positive  and  south 

^  negative) 

(J)  Geocentric  latitude 

U  Argument  of  the  perigee.     The  angle  measured 

eastward  in  the  orbital  plane  from  the  line 
of  nodes  to  the  line  of  apsides. 

l\  Gaussian  determinant 

en  Vernal  Equinox  -  the  point  of  intersection  of 

the  path  of  the  sun  with  the  equator  at  which 
the  sun  is  passing  from  south  latitude  to  north 
latitude  « 

(f)  Angle  at  satellite's  position  measured  between 

-*-  the  line  to  observer5  s  position  and  the  line  to 

the  earth" s  center  (parallactic  angle  at  middle 

observation). 

.,.  Angle  at  observer8 s  position  measured  between 

I  the  line  to  the  earth's  center  and  the  line  to 

position  of  the  satellite  at  the  time  of  the 

middle  observation* 


Q 
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Right  ascension  of  the  ascending  node* 


3.   DETERMINATION  OF  THE  MASS  OF  THE  EARTH 
AND  GAUSS8  GRAVITATIONAL  CONSTANT 

The  existing  relationship  between  t  he  earth  and  the  satellite 
can  be  treated  as  a  modified  two  body  problem  in  which  we  assume  that 
the  two  bodies  are  spheres  and  that  their  mass  distribution  is 
homogeneous*  Under  these  circumstances,  we  can  say  that  they  will 
follow  the  laws  of  celestial  mechanics  and  gravitation.  That  is,  they 
will  attract  each  other  with  a  force  that  is  proportional  to  the 
product  of  their  masses  and  which  varies  inversely  as  the  square  of 
the  distances  between  their  centers  *  Since  the  mass  of  the  satellite 
can  be  considered  to  be  negligible  compared  to  that  of  the  earth  and 
because  of  its  proximity,  it  canbe  stated  that  the  gravitational  attrac- 
tion of  the  earth  is  the  only  attracting  force  influencing  the  orbito 
The  force  acting  between  the  two  bodies  is: 

F  s  k2  %  ms 


where  r  is  the  distance  between  the  two  bodies,  toq  is  the  mass  of  the 
earth,  ms  the  mass  of  the  satellite^  and  k  equals  the  Gaussian  gravita- 
tional constant „     The  numerical  value  of  the  gravitational  constant 
depends  on  the  units  of  mass5  time,  and  distance  chosen©       [2,  Pof>7J 
Kepler's  third  laws 

k2  (me  +  ms)  =  n2a3 
is  a  simple  relation  between  mass,  time,  and  distance  that  may 
conveniently  be  used  in  practical  application  of  elliptic  motion  to 
find  any  one  of  these  three  quantities  when  the  other  two  are  known,, 
The  units  of  mass,  time5  and  distance  can  be  arbitrarily  chosen  and  k 
can  either  be  computed  or  derived  from  observations „     In  applying  the 
units  to  determine  k  for  use  in  observation  of  planetary  orbits  the 
following  units  are  used:  the  sun's  mass  as  the  unit  of  mass,  the 
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ephemeris  day  as  the  unit  of  time,  and  the  semi-major  axis  of  the  earth's 
orbit  (which  approximates  very  closely,  one  astronomical  unit)  as  the 
unit  of  distance o     In  the  case  of  the  artificial  earth  satellite  a 
revision  of  these  units  must  be  made  in  order  to  find  a  realistic  value 
of  the  Gaussian  constant  that  will  be  applicable 0     Bear  in  mind  that 
we  are  theorizing  that  the  earth  is  the  only  body  influencing  the  motion 
of  the  satellite.     This  is  a  sound  theory  but  the  other  assumptions  we 
made  -  that  the  earth  is  a  sphere  and  that  the  mass  distribution  is 
homogeneous-  are  not  absolutely  correct 0     In  fact,  if  the  earth  was 
spherical  and  the  mass  distribution  was  constant  throughout,  then  the 
orbit  of  the  satellite  would  be  nearly  perfect  and  very  simple  to 
predict o     The  plane  of  the  orbit  would  be  fixed,  the  parameters  of  the 
elliptical  orbit  would  remain  constant,  and  the  center  of  the  earth 
would  act  as  one  of  the  focio 

However,  the  earth  is  an  ellipsoid  of  revolution  with  a  bulge  of 
mass  at  the  equator  and  a  non-homogeneous  distribution  of  mass  through- 
out o     These  characteristics  plus  the  effect  of  the  earth's  atmosphere 
and  other  perturbations  have  a  marked  influence  on  the  behavior  of  the 
satellite  in  its  orbital  travels  <> 

The  effect  of  the  attraction  of  the  equatorial  mass  of  the  earth 
is  to  pull  the  satellite  toward  the  equator,  thus  decreasing  the  inclin- 
ation.    Behaving  like  a  gyroscope  the  actual  effect  is  to  cause  the 
plane  of  the  orbit  to  process  westward  for  an  eastward  launch  about  the 
focus.     This  motion  causes  the  perigee  to  advance  to  different  points 
along  the  orbit.     The  perigee  moves  in  the  direction  of  the  satellite 
motion  when  the  inclination  is  less  than  63„°5  and  moves  in  the  opposite 
direction  when  the  inclination  is  greater  than  63  o  £.     Consequently^, 
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the  elliptic  orbit  rotates  about  the  focus  in  the  orbital  plane.     Reduc- 
tions of  the  periods  of  this  perigee  motion  from  the  observations  provide 
data  to  determine  the  shape  of  the  earth*         11,  PjJ 

The  non-homogeneity  of  the  earth's  mass  distribution  causes  per- 
turbations in  the  orbit  of  the  satellite.     By  observing  a  number  of 
satellites  for  long  durations  of  time  we  can  compute  the  deviations 
from  the  normal  gravity  and  thus  obtain  values  for  the  regional  gravity 
anomalies* 

Proceeding  withthe  computation  of  the  Gaussian  constant,  k,  the 
following  units  will  be  used:     the  earth's  mass  as  the  unit  of  mass, 
the  ephemeris  day  (86,^00  seconds)  as  the  unit  of  time,  and  the  semi- 
major  axis  of  the  earth  as  the  unit  of  distance •     Compare  these  with 
the  units  mentioned  earlier  in  this  chapter  that  are  used  in  conjunction 
with  planetary  orbitso 

Note  that  all  computations  are  based  on  the  International  Ellipsoid 
(Hayford,  1910)  which  has  an  equatorial  radius  of  6,378,388.000  meters 
and  a  flattening  of  1/297 °0.     Here  again  rises  an  arbitrary  value* 
There  are  many  reference  ellipsoids  from  which  to  choose  but  the  varia- 
tion in  size  from  one  to  another  is  small*     The  International  Ellipsoid 
has  been  adopted  by  many  countries  and  it  is  the  one  that  is  used  by 
Smithsonian  Astrophysical  Observatory  in  their  computations,  so  I 
decided  to  use  this  as  the  reference  body. 

The  international  gravity  formula  adopted  by  the  International 
Union  of  Geodesy  and  Geophysics  in  1930  is     p*  P«7M    s 

7s    Z.  <1  +  J3   si*2<£   +     £    sin220  ) 
where      [ e  »  gravity  value  at  the  equator,     (f)  s  latitude  of  station* 
Q   z  gravitational  flattening,    €  *  theoretical  constanto    The 
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coefficients  of  the  formula  are: 

7s  978.02*9(1  +  0.0Q5288U  sin20    -  O.OOO0G$9  sin22(ji>)c  jn/sec2 

The  gravity  value  at  the  equator  is  the  resultant  of  the  gravitational 

attraction  and  the  centrifugal  force  due  to  the  earth's  rotation. 

There  fore , 

Te  -  w2a  ■  978 ,01*90  cm  /sec2 

where  Te  =  attraction  of  the  earth  at  the  equator  and  w2a  =  centrifugal 

force. 

If  w2a  s  (-=-)2rjWhere  w  is  the  angular  velocity  of  the  earth, 

T  is  the  sidereal  day  in  seconds. and  r  the  radius  of  the  ellipsoid  them 

(4j2  )2r  -  3#38o#i5  cm/see2 

Then 

Te  -  3»3805  =  978.01*90 

Te  =  978.01*90  +    3-380^  a  98l.l*29S  gal 

where  a  gal  (from  Galileo)  equals  a  unit  of  acceleration  equal  to 

1  cm/sec2 • 

The  mean  density  of  the  earth  is  computed  from  the  formula    I  5>  1  s 

where  P©  ■  Te/a  =  attraction  of  earth  at  the  equator  divided  by  the  semi- 
major  axis,     G  s  Newton's  universal  constant  of  gravitation  which  equals 

ft  I  ~s 

6.673  i   0.003  x  HP     cgs,  e     =  the  secondary  eccentricity,  and      Q,^  s 
the  mean  density  of  the  earth.     To  solve  the  equation,  two  of  the  three 
quantities  Te,  a,  or   Urn  must  be  known.     In  our  example,  the  Te  has 
been  determined  and  the  value  of  the  semi-major  axis  of  the  Interna- 
tional Ellipsoid  is  used*.     The  computed  mean  density  is  (see  Table  1): 
CTffi  s  £.£13001*37  gm/cm3 
The  volume  of  the  International  Ellipsoid  is  computed  from  the 
formulas 
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V  =  k/3    Tf     a2b  =  1.083320781  x  1027  cm3 
where  t>  =  6,356,  911  «9li62  meters,  the  semi-minor  axis  of  the  ellipsoido 
Using  these  values,   the  computed  mass  of  the  earth  is: 
/v/    =  \/Cf    =(1.083320781  x  1027)  (5.51300207)=  5.97235220  x  1027gm« 
To  obtain  the  proper  unit  for  k  in  terms  of  mass  of  the  earth  for 
unit  mass,  time  in  days,  and  distance  in  terms  of  the  earth's  semi- 
major  axis  the  following  conversion  is  made: 


■■?'■'     cm3  = 
see* 


cnP 


\2   (c  aio'it.oo  -»-  in27' 


gm  sec 
and  k  =  107.0731376  and  k2  =  llU61i.6568 


gm=     n  JSh^JZ  (5.9723522  x  ,10*  ,)  =  k 
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TABLE  1 
COMPUTATION  OF  VOLUME,  DENSITY,   AND  MASS  OF  THE  EARTH;   COMPUTATION  OF  k 
DIMENSIONS  OF  THE  INTERNATIONAL  ELLIPSOID 


a    a  1/297 

a      -  6.378,388.0000  meters 

b      =  6,356,9H.9U62  meters 

e      =  0.0819918898 

e2    =  O.OO6722670O 

e      =  0.0822688896 

/-?  »  0.0067681702 

e,;J  =  0.0005568098 


li>  2  -  1«0Q67681702  =  i808.10066 
.0005$6BO9» 


75" 


ARCTAN  e  =  e  - 


/-J?     e73       e'S-e'7 


3    ^~T     "7 
...  =    0.082081|036$ 


+ 


J: --. *  0.0817158230 


(-^L)2r  =  (2  ^  3.1la592653^2  6,37.8,38800.  =  3-3805U56UOO 
7"  d6,ii00 

V  -  U/3  7f  a2b  =  1.083320781  x  IO27  cm3 


tf»  = 


■27re 


1+e' 


7TT 


^3 


arctan  V  — 


5.51300U37  gm/cm3 


l+e' 


n 


(a) 


M  =   VCm=  5-972352200  x  IO27  gm 


k  =  6-L,  «  -  G 
sec* 


cm- 


gm  sec 


gm=   G  (86^00)2(5,972352  x  1027 
(6,37&,38800p 


107.0731376 


U.     DETERMINATION  OF  THE  ORBIT 

Since  the  orbit  of  the  satellite  is  defined  by  six  parameters  then 
six  independent  quantities  must  be  given  by  the  observations  in  order 
that  the  parameters  may  be  determined.     A  single  observation  gives  two 
quantities,  the  angular  coordinates  of  the  body.     Therefore,  three 
complete  observations  are  sufficient  to  define  the  orbit.        U,P»19lJ 

The  apparent  positions  of  the  observed  satellite  are  obtained 
by  measuring  its  angular  distances  and  directions  from  a  background 
of  stars  on  the  film  of  the  Baker-Nunn  camera.     The  stars  on  the  film 
can  be  identified  from  catalogues  containing  their  right  ascension  and 
declination.     With  these  coordinates  known,  then  the  right  ascension 
and  declination  of  the  satellite  cante  computed.     Catalogues  containing 
the  observed  values  of  time,  right  ascension,  and  declination  reduced 
from  the  camera  film  are  published  about  six  months  after  the  observa- 
tions were  made.     These  reduced  values  are  quite  accurate.     However,  a 
more  precise  reduction  of  the  films  is  made  by  the  Smithsonian  Institu- 
tion Astrophysical  Observatory  Photoreduction  Division.     These  results 
are  published  and  made  available  about  18  months  after  the  observations 
are  taken* 

In  this  paper  the  computation  of  the  orbit  of  Explorer  I  (19E>8  Alpha) 

satellite  will  be  reduced  from  observations  of  right  ascension  and 

declination  from  the  Baker-Nunn  Camera  Tracking  Station  at  Olifants- 

fontein,  South  Africa.    Date  of  the  observations  was  18  December  I960* 

The  observed  data     I  7,  p.lOj  s 

Time  (UT)  C£  0 

22*1    Vf    10s. U*  Sh  07m  06S  -1*7°  22 '  00" 

22h    31m    0U°.li2  7*J  28m  &s  -£j.°06«Q0H 

22h    32m    03SJL0  8ft  33     00s  ~k9°  2$'00n 

These  data  are  not  those  precisely  reduced  by  the  Photo-reduction 

20 


Division  but  the  precision  of  the  observations  is  as  follows: 
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7,  Po3 


Standard  Error  in  timing  (    Cft    )  .0Q£S     ■*     tffc     =      .02s 

Standard  Error  in  direction  (    (f^    )         2f.7      -*     V*      ^      3'.£ 

The  geographical  coordinates  and  height  of  the  tracking  station  are 

[8,  p.7]  s 

0  =  25°  57"  3U".  70  S 
A=  28°  1U1  51"  .10   E 
H   =    lf&U  meters 

There  are  two  classical  methods  available  for  determining  the  orbit- 
one  developed  by  Laplace  in  1780  and  the  other  by  Gauss  in  l801o  Both 
methods  were  devised  for  use  in  determining  the  orbits  of  comets  and 
planets  about  the  sun  but  are  also  applicable  for  use  with  artificial 
earth  satellites.  Since  the  development  of  these  two  methods  many 
papers  have  been  written  on  the  theory  of  orbital  determination  but 
very  little  that  is  theoretically  important  or  radically  different 
has  been  added  to  the  methods  of  these  two  men.  The  Laplacian  method 
is  not  as  complicated  as  the  one  devised  by  Gauss  but  it  has  two  weak- 
nesses. It  is  inaccurate  when  the  intervals  between  observations  are 
short  and  it  assumes  that  the  observer  is  in  orbit  about  the  focus  of 
the  ellipse*  The  Gaussian  method,  even  though  it  is  an  approximate 
method  was  chosen  to  determine  the  orbito  There  are  other  more  modern 
solutions  available  but  these  are  not  as  accurate. 

The  first  step  in  the  calculations  is  the  conversion  of  the 
Universal  time  of  the  observations  to  sidereal  time.  This  is  required 
in  order  to  introduce  the  vernal  equinox  into  the  coordinate  system* 
The  sidereal  time  at  any  instant  is  the  hour  angle  of  the  vernal 
equinox • 

The  origin  of  the  coordinate  system  is  at  the  center  of  mass  of 
the  ellipsoid.  The  x  -axis  points  toward  the  vernal  equinox,  the  y-axis 


22 

is  directed  toward  a  point  90°  eastward  from  the  vernal  equinox,  and 
the  z-axis  coincides  with  the  north  pole  of  rotation.     The  x-y  plane 
lies  in  the  plane  of  the  equator. 

The  rectangular  space  coordinates  of  the  observer  are  computed 
from  "the  given  latitude  and  longitude  of  the  tracking  station.     Then, 
by  changing  the  signs  of  these  values  the  coordinates  of  the  center  of 
the  earth  are  obtained.     The  observer's  position  is  computed  as  follows: 

xj  =  (N  +  H)  cos  0    cos     Q  3.  (j  =  1,  2,  3) 

yj  =  (N  +  H)  cos  0    sin    0  lj 

zj  =  (N  +  H)  sin  0     (1  -  e2) 
The  position  of  the  observer,  the  earth's  center,  and  eventually  the 
position  of  the  satellite  are  all  converted  to  rectangular  space 
coordinates  to  provide  a  more  convenient  and  a  common  system  of 
coordinates. 

The  next  step  is  to  compute  the  direction  cosines  from  the 
observer's  position.     These  are  functions  of  the  observed  (apparent) 
right  ascension  and  declination  and  are  obtained  from: 

f,  =  cos  0,  cos  CX  .  U  =  1,  2,  3) 

gj  =  cos  0  j  sin  CX  * 

h-j  =  sin  6  j 
The  orientation  of  these  axes:   f-axis  points  toward  the  vernal  equinox, 
g-axis  perpendicular  to  f-axis,   and  the'  h-axis  parallel  to  the  axis  of 
rotation  of  the  earth  and  perpendicular  to  the  f-  and  g-axes.     The  direc- 
tion cosines  give  us  a  line  in  space  from  the  observer1  $  position  to  the 
satellite.     The  computations  of  these  three  steps  for  the  three  observa- 
tions are  contained  in  Table  2» 


TABLE  2 
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COMPUTATION  OF  SIDEREAL  TIME,  COORDINATES  OF  EARTH  CENTER,  DIRECTION  COSINES 


Date 

Sidereal  Tine  At  V1  DT 

UT  Time  Interval 

UT  To  Sidereal  Time 

Greenwich  Sidereal  Time 

A  (Add  East,  Subtract  West)   +  1,-52-59.1+07 
Local  Sidereal  Time  (19  Dec)       6^12^163599 


18  December  I960 
5-1+6-25.1+19 
22-29-lO.liiO 
3-1+1.633 

28-19-17. 1$2 


son. 


e2ski20 

1  -  e2sin2j 

(l-*e2sin^ 

a/(l-e2° 


Bin 


H 

H/a 

N/a 

N/a  +  H/a 


N+H 


H)  cos  (£>  cos  £?£ 
H)cos  (7)  sin  C/| 
H)sin  (?)  (1-e2) 


-  x 

3  7 

s  z 


93°0l+'08".985 

.99856561+ 
-.05351+139 

^_o7m_o6s 
76°U6'30» 
.9731+7917 
.22877565 

it7°22'00»  S 
-.73570317 
.677301+10 


cosQC      =  f 

slaQC'     =  g 
=  h 


+  .01+818189 
-.89860916 
+  .1+3518070 

+  .151+95069 
+  .65931+H+3 
-•73570317 


18  December  I960 

5-1+6-25.1+19 

22-31-OU.1+20 

3-1+1.91+6 

2^-21-11.78$ 

+  lr52-59.1+07 

93°32'li7M.68o 

.998081+78 
-.06186091 

7h_28m-51is 

112°13,30» 

.92570563 
-•3782UJU7U 

5l°o6«oo»  s 

-.77821+315 
.62796306 

25°57»3U».7  S 
-•1+3773790 
.1916H0+7 
.89910263 

6,378,388.0000 
.00672267 
.99327733 
.00128816 


18  Dec.  I960 

5-1+6-25.1+19 

22-32-03.100 

3-1+2.107 

28-22-10.626 

+  lr  52-59.  1+07 

mgomm 

93°l+7'30".l+95 

.99781091+ 
-.06613118 

8h_33m^)os 
128°15«00M 
.78531693 
-.61909395 

l49°25'00M  S 
-.7591+6058 
.65055333 


.99871181+  ■*■ 

.99935571  -*■ 

6,382,500.18  meters  * 


151+1+  meters  **- 

.00021+207  *~ 

1.000614*71  <=- 

I.OOO88678  <" 


+  .05566863 
-.898176U3 
+  .1+3518070 

-.23752372 
+.5813089U 
-.77821+315 


+  .0595111+1+ 

-.89793000 

+  .1+3518070 

-.1+0275363 
+.51089051+ 
-.7591+6058 


2k 


The  next  step  is  to  determine  the  distance  of  the  satellite  from  the 
observers  position  and  from  the  center  of  the  earth.     This  is  done  by 
solving  the  triangle  formed  at  the  time  of  the  middle  observation  between 
the  center  of  the  earth,  the  observer,  and  the  satellite.     To  solve  the 
triangle  the  satellite  is  "stopped"  in  its  orbit  at  the  time  of  the 
second  or  middle  observation,     (T2)« 

Equator 


Figure  k» 

The  geometric  situation  is  depicted  in  Figure  k»     Since  SOE  form  a 
triangle  at  T2  the  quantities  r2  and  R2  must  satisfy  the  equation 
[U,p.20l]    : 

r|  -  r2  +    p  2  -  2rp      cos  ^ 
This  equation  treats  the  geometric  situation.  In  the  Gaussian  method 
a  formula  relating  the  two  unknown  quantities  r  and  R  applies  dynamically 

^P,235]    S  -     "      *  ,  t2Kl 


r2 


~T*T 
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where  K,  LA  t  and  Ki  are  determinants  involving  the  coordinates  of  the 
earth's  center  and  the  direction  cosines  from  the  observer's  position. 
These  determinants  and  their  values  are: 


A 


K 


ICi-  - 


"  fl  f2  f3~ 
gl  g£  g3 
hi  h2_  I13 

3 

-  0.0020U368 

*1      xi  +  x3  • 

gi    yi  +  73 

_  hi      zi  +  Z3 

-  2x2         f  3  ' 

-  2y2        g3 

-  2z2        h3 

"  f        xi  +  x3 

gl      7i  +  73 
hi      zl  +  z3 

h3 
h3_ 

=    +  0.00018897 


=      -0.li57l6U31 


The  solutions  of  these  determinants  are  shown  in  Table  3o 

The  two  unknowns  r2  and  R2  are  expressible  in  terms  of  a  single 
unknown,  the  angle    4>    ,  by  means  of  the  equations 


P2  cos     ¥   2  =  x2  f2    +     72  g2   +   z2  k2 


Nsin  m 
[\|  cos  m 


>  = 
=  P 
*  P2  c°s 


U.p.235]   ! 
(   0  -  Y  =  IT) 


¥..• 


2  san  1  2  =     +  .U8507360 


A 


The  computed  values  for     ip,  N,  and  m  are: 

M/2    =    IJO0  58'  12«.06H 
N      =      I  0.91981i731 
m      =    lt#°  10'  26" .395 


The  solutions  for  these  quantities  are  contained  in  Table  k* 
N,  m,  and  M  are  intermediate  parameters  used  in  the  solution  for  the 
angle    $    ,  where     \k$  p.212]    : 

M=    k     ^    {%   N,  |j^,4  0.163U1&3 
The  sign  of  N  is  chosen  so  that  M  shall  be  positive. 
The  equation  for  the  solution  of  the  angle     4?  is  |U,  P.213J  : 
sin  *  Q    =    M  sin  (  (J)  +  m) 
which  is  a  transcendental  equation.     The  quantities  M  and  m  are  known  and 


M  is  positive.     One  solution  of  the  problem  results  in  the  position  of 


TABLE  3  26 

SOLUTION  OF  DETERMINANTS 


A                                  ?1  82  h3  =  -  «068  ^°7  82 

A                                f£  g3  hi  -  .089  276  57 

^3  gl  h2  a  *206  661*  1^ 

f3  g2  hi  -  .172  2U5  98 

fl  g3  h2  =  -  .061  607  9k 

f2  gi  h3  -  .118  938  5U 


K        fi  (yi  +  73-2y2)h3  =         .000  021  92 

(X1+X3  -  2JC2)g3hi  =         .001  369  62 


f3gi(ai+  z3-2z2 )  =  0 

5  Tyi+y3-2y2)  hi        =     --000  °- 

fl  g3(zi+Z3-2z2)  s  0 


Tyi+y3-2y2)  hi        =     --000  °5$  20 

S3(zi+Z3-2z2)  s  0 

(xi  +  X3-2X2)  glh3  3         .001  82U  68 


Ki     fi  (21+73)  h3  =  .211  iau  83 

(xi+x3)  g3  hi  =  -.OliO  U78  02 

*3  SI  tzi  +  z3)  ■  -.231  126  3k 

*3  (yi+  y3)  hi  s  -.532  327  35 

fl  S3  (zl+z3)  =  'O68  900  28 

(xi+x3)  gi  h3  =  -.053  926  77 


x2  f  2 

y2  s2  = 

z2  h2  — 

x2  f2  +  ?2  22+  z2  n2  = 


TABLE  1* 

SOLUTION  OFx/2*  N,  m. 

-.013  222  62 
-.522  117  99 
-.338  676  1*0 
-.871*  017  01 
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P. 


cos 


% 


(x*  +  y?  +  zg)2 


.003  099  00 

.806  720  90 
.189  382  2l* 
.999  202  11* 
.999  600  99 


*2f2  +  y2g2  +  z2h2/P2=   -  -87U  365  89 
sin 
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P2  sin    Y2 

P2  cos    y2  -  JL 

N2  sin2  m 
N2  cos2  m 
N2(sin2m-f  cos2m) 
N 

P2  sin  X  2  / N 


150°  58'12M.061* 
.1*85  267  23 


P2  cos 

A 


A' 


P\ 


N 

U  A  P  I  N  si^3^ 

t2  Ki 

U  API  Nsin; 


-1.  3M/ 


Solution  for  t: 


3 


t  = 


k(T3     - 


Ti 

Tl 
k 

Tl) 
t 


=  .1*85  073  60 

-.781  551  U6 

=  .235  296  1*0 

=  .610  822  68 

=  .81*6  119  oa 

=  .919  81*7  31 

.527  3U  1*3 
-.81*9  653  U7 

=      11*8°  10«26».395 

=        -.002  01*3  68 
=  .998  803  U5 

.111*  272  81 
.919  81*7  31 
*  .011  1*85  88 

=        -.161  Uk  31 
=         -.000  858  2f* 

=        -.005  250  93 
.163  W6  33 


k  (Tq  -  Tl) 


T3  =      22h  32™ 

T-i  a      22h  2^ 


03s. 10 
lQS.H* 
2111    52s. % 
.002  001  85  days 
107.  073  1376 
0.211*  3hh  36 
0.107  172  18 


P. 


2 

cos     X ; 


N  sin  M 
N  cos  M 


sin  m 
cos  m 


M 


=       2t 
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the  observer  and  is  not  applicable.     That  is,  if  r  =  0,    P  =  R.     The  proper 
solution  for  the  angle  Q?    must  satisfy  the  inequality     [ti,  p. 213     : 

In  our  problem,    X  equals  l£00  58,12",06U  so    :r    must  be  less  than 
2°°02'.     The  solution  maybe  obtained  graphically  be  plotting  the  curves 
U,  P. 213]    : 

y2  -  M  sin  (  (T)  +  m) 
Once  the  approximate  value  is  found  from  the  intersection  of  the  two 
curves  then  the  exact  solution  may  be  found  by  numerical  trials  or  by  a 
differential  formula.     The   solution  for  the  angle  is  lengthy  and  time  - 
consuming.     The  approximate  value  from  the  intersection  of  the  curves  on 
the  graph  was  found  to  be  approximately  23  °»     The  exact  value,  solved  for 
by  numerical  trials  was: 

(J)    Z    23°  l£»  U9'M6 
Uith  this  value  known  we  can  now  solve  for  the  two  unknown  sides  of  the 
triangle  SGE  in  Figure  h9  which  will  yield  the  distance  of  the  satellite 
from  the  observer's  position  and  the  distance  of  the  satellite  from  the 
earth's  center  for  the  second  observation.     The  formulae  used  and  the 
results  are     I  kf  P«  25&1    s 
R2 


=    Pgain^ 


sin(F)      T 

Pgsinj^ 


«    1.228  1^0  03 
)-    0.2&  280  87 


sin 


The  observer  -  satellite  distances  for  the  first  and  third 
observations  are  given  by      k9  P»  25£J : 


glS3 


ri  = 


xXf3 
yiE3 


x2f3 
7263. 


♦*$ 


x3f3 
7323 


+  r2 


m 


f2f3- 

g£g3 
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V3 

&Lg3 

where 


=BteW 


flx2 

S172 


+ 


fxx3 


+  *Z 


w 


h*2 


0^2  =  1-3  + 
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i^N+z 


t8 


=    1.9U3  20695 


»ai*yy 


^i*1 


...  =      0.5U*6l323 


(}>3)  =  1  -      1.50993233 

2      2t  +  1*1^ 


where 


t=  k  (^  -  sj  = 


0.10717218 


q  =    k  fT1  f  |3J     =      -0.031*1*5185 

The  time  of  the  second  observation,  T2,  is  the  origin  of  time.     Upon 
substituting  these  values  into  the  above  formulae  and  solving  the 
determinants,  the  observer  -  satellite  distances  for  the  first  and 
third  observations  are  found  to  be: 

ri  =  0.23531756 

r3  =  0.27239982 
The  detailed  solutions  for  R2>  *1»  T2>  r3  a*1**  Q  a3*6  listed  in  Table  5» 
The  next  step  is  to  determine  the  geocentric  space  coordinates  of 
the  satellite  for  the  time  of  each  observation  as  follows     ll*,  p. 256  : 

r  0.23531756  0.251*28087  0.27239982 

X  =  r  f  -  x  -0.01171927         -O.II606637  -0.1692211*5 

Y  =  r  g  -  y  + 1.05376378       +  1.0U599217         + 1.0370961*9 


TABLE  5 
DETERMINATION  OF  R2,  r2,  q*  ri,  Yy 
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# 


^2  sinx*    2 

H2  sin      x  2/sin 

I7  = 

<£  +   $>*  = 

sin  f  ^  2   -f  <±)  )  _ 

p2  sin  (yF+flb'L  - 


T3         »  2^  32m  03s  .10 
T2         ■  2£&  31m  0lis.U2 

T3  -  T2   =  (+)   58s.68 


0.1*85  073  60 
0.39U  962  82 
1.228  ISO  03    ■  R2 

150°  58' 12". 06 

17ltu  Hi'  01"  o 22 
0.1COU7158 
O.lOOltflltf 
.25U28087  =    r2 

T2  *  22**  31m  0US.U2 
Ti  =  22h  29m  10s. Hi 
m 


T2  -  Ti 
T3  -  T2 
Ti  +  To 
T1  +  T3 


•  (-) 
■  (+) 
»(-) 
(days) 
k 
k  (T1+T3) 
k  (T1+T3)/2 


T2  -»j  =  (-)     lm  51is.28 

lm  51ts.28 
58s. 68 

^s^5 

=    -0.00061*352 

=   107.0731376 

a  -Oo 06890371 

=    -0.031445185  =    q 


t 

2t 

t2 

R 

R2 

2R2: 


0.10717218 
0.21i43Ui36 
0.0111*8588 
1.22815003 
1.8521*8316 
3.70it96632 


hR2 

q 

q/2t 

o  <*/** 

t2/2R2^ 

t2/^23 


7.1*0993261* 

-0.031*1*5185 

-0.16073131 

-0.3211*6262 

0.00310013 

0.00155007 


1  -  q/t  +  t2/2R23  =      1.321*56275      1+q/t  +  t2/2R23  =      0.68163751 
l/|+|  +  |L      =     iLU-  2.J31O1079 


1  - 


I + W/1 + 1 + £?  =  M  =  1,91°20695 


+l^/1-^+^ 
/*** 


$8  = 


0.511461323 


+  T3ET 


=     1.50993223 


TABLE  5  (Continued) 
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X1  53 

n  % 

xig3  -  n^3    = 


x263 
y2*3 
x2g3  -    72£3    ' 


T2757    [72^_ 


yi63 


x2f3 

y2&3 


c3|3 


X-: 

*3g3  -  y3^3 


x3£3 
7383 


(1*2) 

T57f) 

f2?3 

82*3 

f2g3-  g2f3 


%$} 


81*3 
^183-  81*3 


flf3 
8183 


[x3f3 

73^3 

f2f3 

g2*3 

\£j£z 

g283 

&f^ 

g!83 

0.02U61567 
0.36191810 

-  o.mwm 

0.02814*058 
0 .361711382 


-  0.97791*863 

0.030^0383 
0»36l614*57 

-  0.3312U07U 

-  0.6U366931 

-  0.12131*862 

-  0.231*121*29 
0.11277567 

0.08Umo60 

0.07916281* 

-  0.26555215 
0.3*4*711*99 

0.0811171+9 


rx      =    .08lll7l*9/.3l4*7U*99    =    0.23531756 


Sl*l 
*i3tt  -  81*1 


'Wl 


fixi 


TT727 

gl*2 
f 172  -  gl*2 

(1,3) 


flxl 
8121 


fix2  1 
gl72 


flx2 

8172 


-  0.13 921*011 
0.03176832 
0.1710081*3 


-  0.08800320 

-  0.13917306 
0,03670li63 

-  0.175B7755 

-  0.26556339 


TABLE  5  (Continued) 
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*U3 

21x3 
£373  -  gix3 

3 

"fix3" 
2173 

flf2 
21^2 
fl22  "  21^2 

= 

"flf2" 
glg2 

r2  Ws 

"*lV 

2122 

"flf3" 

2123 

r3 

-  0.13913ii87 
0.03923836 

-  0.17637323 


0. 090071*22 

-  0.15660923 

0.2l*6683US 

0.091*71331* 
0.09390030 


r3      -       .09390030/.3l*l*711i99    =      0.27239982 
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2  =rh-z  -  0.60830i£7  -  0.63307305        -  0.6U205763 

R    =V(x2+Y2+22)  1.2l679k93  lo228l£003  1.231U39U0 

There  is  a  check  for  accuracy  of  the  computations  up  to  this  point 
in  this  last  step»    Note  the  agreement  in  the  value  for  R2  computed  from 
the  space  coordinates  and  that  determined  from  the  formula,  {P2  s^21    ±2)/ 
sin  4r  ,  in  Table  £• 

The  inclination  of  the  orbital  plane  and  the  intersection  of  the 
orbital  plane  with  the  plane  of  the  equator  (XY  plane)  can  now  be  solved 
using  the  geocentric  space  coordinates  in  a  ratio  of  determinants ->     This 
solution  will  produce  the  first  two  parameters  of  the  orbit  that  we  are 
seeking:  i  and  \  I   \k9  p. 2£6     : 
A:  Bs  C     = 


Y222 


Z1X1 
Z2X2 


XI  Yi 
X2  Y2 


See  Table  6  for  the  evaluation  of  the  ratio  of  these  determinants  and 
the  detailed  solutions  f or  i  and  i  :  • 
Then, 

cos  i  ■  C_ 


sin 


cos 


Q 


Q 


Va2  + 

B2  +  C2 

sin  i 

= 

A 

Va2 

+  B2  +  C2 

sin  i 

- 

~B 

i  =  32°  3U»20«.2l£ 


v 


Q  =  206°  00«27Mo603 


'A2  +  B2  +  C2 

If  the  inclination  was  greater  than  90°  the  sign  of  the  previous 
two  equations  would  be  reversed.  The  values  i  and  ^  I  uniquely  determine 
the  position  of  the  orbital  plane  in  relation  to  the  equatorial  plane 


of  the  earth 


U,  p. 11*6 


The  next  step  is  to  find  the  auxiliary  parameter  u  for  each  observa- 
tion. Solution  of  this  value  is  necessary  in  order  to  find  the  "longitude" 


TABLE  6  3U 

EVALUATION  OF  DETERMINANTS     A:  B:  C  5  SOLUTION  FOR  i  AND  W  . 


*L*2 

Y2Z1 
*lZ2  -  Y2ZX 

zlx2 
Z2Xi 

Z1X2  -  Z2Xi 


X!T2 
X2Y! 
XiY2  -  x2n 


IlZl 
Y2Z2 

Zl*L 
Z2X2 


X2Y2 


A2  +  B2  +  C2 


C/^\/A2  +  B2  +  C 


-0.667109l;5 
-0.63628182 
-0.03082763 

0.07060370 
O.OQ7U1915 
0.06318U55 

-0.01225826 

-0.1223065k 

0.1100U828 


0.13058811 
0.8U271286 


=      B 


-       cos  i 


sin  i 
A2  +  P  +  Cd 


A_ 


v 


A2  +  B2+  C2/ sin  i 

/A2+B2  +  C2 
B 


V 


A2  +  B2  +  C2/sin  i° 


ft    = 


32°3V20".2H5 

0,53836329 
-0.23606766 


-0.U38U9138 
O.U838U632 

-0.89873535 
2O6°00»27M.603 


-   sin 


Q 


COS 


ft 


3* 


of  the  satellite  from  the  position  of  nodal  passage.  The  longitude  of 
the  satellite  from  the  node  is  called  the  argument  of  the  latitude • 

S3 


Figure  £. 
Refer  to  Figure  3>»  The  formulae  that  define  ui,  U2,  u^  are  U,p,238 

cos  Ui  =  cos  b.  cos  (I*  -  W  )        (j  ■  1,  2,  3) 

».  «       j       j      j 

sin  uj  cosi  s  cos  bj  sin  (lj  -  W  ) 

sin  u.'  =  sin  b.»  /sinj_ 
To  determine  the  parameters  a,  e9  and  w;we  must  solve  for  the  ratios 
of  the  triangles  formed  at  the  time  of  each  observation.    With  the  geo- 
centric coordinates  and  the  range  of  the  satellite  from  the  earth* s 
center  known  the  above  equations  may  be  substituted  with  the  following 


expressions 


k,  V.2$3S2$1 


sin  i  sin  u*  - 


cos  x  sin 


UJ 


cos  u- 


Then  for  abbreviations 


cos  i  I 
sinO    + 


RJ 


Rj 


j  -  (1,  2,  3) 
sin  i  i 
cos  j  t 


(f    »  u3  -  u2    =  2°oiiU97583  s  0.0li2620$0  radians 
(f2  a  u3  -  ui    -  7°.U^16083  -  0.13011710  radians 
CT3  .  u2  -  ui    a  £°.0131&S0Q  s  0.087U9658  radians 
Then  the  ellipse  parameter  p  is  determined  and  used  in  place  of  the  semi- 
major  axis,  a. 
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The  equation  that  defines  p  is     h»  p. 2^7     : 

kV^1  (t3-.tx)  s  r|  cf23  +  r^2  (f2(2cr3~gr1)+R3cr2(2g;~  C/3) 

6  01    (J 3  6    ^3  6    cT1 

Using  this  formula  the  value  for  p  was  found  to  be: 

p  =    0.830&3& 
With  the  values  of  p  and  u  now  known  the  eccentricity,  e,  and  the 
argument  of  the  perigee,  w,  can  be  found  from    U,  p. 2^7   : 

e  sin  (ul  -w)  =  i^p-Ri)  cos   V*  -  R-i    (pgR-Q 

Rl  R3  sin    (j  2~ 

e  cos   (u^-w)     a  p  -  Ri 

Then  the  semi-major  axis  of  the  orbital  ellipse  is  defined  by 

a  =     (1  -  e*) 

The  computed  values  for  these  three  orbital  elements  are: 

e  =  0.3281^828 
a  =  0.93078812 
w  =     82°^,U2».337 

The  detailed  computation  of  these  parameters  and  the  auxiliary  para- 
meters ux,  U2,  U3,   and  p  are  listed  in  Tables  7 »  8,  and  9» 

The  next  and  final  orbital  parameter  to  be  computed  is  the  time  of 

epoch,  T0  (time  of  perigee  passage )o     This  value  is  computed  from  Kepler's 

equation: 

M  =  E  -  e  sin  E 

after  the  true  and  eccentric  anomalies  are  determined  for  the  three 

observations.     Values  for  To  are? 

Observation  #1  Observation  #2  Observation  #3 

December  l8«,9m5H9U2  18.911*53182  18.91U5S056 

Detailed  computations  of  T0  are  listed  in  Table  10© 

To  verify  the  results  of  the  computations,  the  orbital  parameters 
were  determined  again  using  data  published  in  Smithsonian  Institution 
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Astrophysical  °bservatory  Special  Report  #78  dated  25  October  1961,     This 

publication  contains  both  the  mean  elements  and,  the  smooth  elements  of 

Explorer  I  (1958  Alpha)  satellite  based  on  82  observations  made  from  the 

various  Baker-Nunn  Tracking  Stations  for  "the  period  16  December  through 

29  December  I960.     The  published  smooth  elements  ares 

To  a  37291.0    Modified  Julian  Days  (Smithsonian  Days) 

w    =  (U*0°.6  ±2)+7°.275  (t-To)  +  0°.000U27(t-To)2+0°.22U  cos  w 

Q   =  (186°. 278  ±9)-Uo.896(t-To)«0o.000298(twTo)2  +  0°. 00887  cos  w 

i    =  (33°. 210  ±3)-0°.000^U(t-To)+0°.50xl0-,^(t-T0)2  -0°.OOii73  sin  w 

e     =  (.09971   ±8)  -0. 000202  (t-T0)  +  0.30  xlQ-7(t-T0)2+  0. OOOU5519  sin  w 

M    =  (.9U016  ±8)     (13.U3S13U  ±5)  (t-T0)  +  (0.0002839  ±6)  (t-T0)2 

-(.11  ±1)  xl0-5(t-To)3  -< .000631  xlO-3)  cos  w 

Standard  error  of  one  observation:     m  =     ±5°<»U3<> 
The  angles  are  expressed  in  degree s,  e  in  terms  of  the  earth's  semi-major 
axis  (6,378,388  m.),  and  M  in  revolutions.     The  Smithsonian  Day  for 
December  0,  I960  is  37268.0.     Using  the  time  of  epoch,  T0  «  December 
18.911*53182  from  the  middle  observation  of  our  calculations  and  redesigna- 
ting it  as  t,  the  term  (t  -  T0)  will  have  the  following  value: 

Date  Smithsonian  Day 

December  0,  i960  37268.0 
Date  of  Observation  +  18.911*53182 

Date  of  Observation  (t)  37286. 911153182 

Epoch  (T0)  37291.0 

(t  -T0)  -li.68&6818  days 

Detailed  computations  of  the  elements  are  contained  in  Table  11.     The 

computed  and  published  parameters  are: 

Published  Computed 

i  33°13«OO.ii02  32o3U'20.2U5 

Q  206°l6«il3".067  206°00«27Mo603 

e  0.10011028                               0.32815828 

a  1.17098895                               0.93078812 

w  110°57»53M.75U  82°5U'U2M.337 


TABLE     7 
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sin  I 


(Zj/Rj)/sin  i  =  sin  u 


COMPUTATION  OF  THE  AUXILIARY  PARAMETER,  u 

Observation  1     Observation  2  Observation  3 

-O.U9992366   -0.515U6882   -0.52138792 
0.^3836329    Oo^3836329    0.53836329 
-Oo  92859909   -0.957V739U   -0.968U6856 


cos  -1  L 

Y/R  sinr) 

X/R  cos  * L 

Y/Rsin  Q  +     X/R 


u3""u2  =  ^1 
Ol  (radians) 

U3-U2  r  0 2 
($2  (radians) 

u2-u1  =  $3 
03  (radians) 


cos 


ft 


0. 86601592 

-0«00963126 


cos  u 


°0.3797U052 
0. 00865595 
-0.37108U57 


0.85168110 
-0.09U5050U 
-O.U38U9138 
-0 o89873535 
-0.373U5U82 

0.081i93502 
-0.28851980 


0.81+218232 
-0.137U1760 


-Oo36928969 

0.12350205 

-0.2U57876U 


2it8°13°02M.81i5    253°l3'50w.3ll  255°Uo«21".U2U 


2026«31".113 
0.0U262050 


7°27fl18».579 
0.13O11710 


5°00'U7"«U66 
0.087U9658 


TABLE    8  39 

DETERMINATION  OF  AUXILIARY  PARAMETER,  p. 


kVF1     (T3-Ti)  *    Eg2     CT23         Rl2  CT2(2  0^-  (T-, )  ,  R,2  Cf»(2   Cfe  0Q 

60i  CT3  6  cT3  6    (^ 


Rl2  1.U8058990 

R22  1.508352^0 

R32  lo5l6Ut300 

(M  0„0022029U 

fl/fet.  Oo0033228l 


0.022371*90 

R22     09V(6  CTi  CTo)  Oolii85o6m 

R12    rf?  0.19265006 

(2  <£i  =  0,)  0.13237266 

6  o^3          t       ^          ^  Oo52U979U8 

Rl2  0?(2  03-  tfi)/6  03  O.OU857637 

R32  fti  Ool973l5l7 

(2(Ll-  03)  -0.00225558 

6„(X                 _^           .  0,25572300 

Rj  02(2  dV-  JO/6  Oi  -o.ooi7l*oUo 

k  Vp-  (T3-Ti)s  #+»+*    s  Ool953U211 

K  (T3-T1)  0o  211*31*1*36 

*+»+*A(Tq-T!)  ■  y^~    a  0.9113U710 

P       J  0.83055352* 


TABLE  9  1*0 

COMPUTATION  OF  THE  ORBITAL  PARAMETERS  es   a,  w, 

e  sin  (ui=w)  ■  R^  (p-Ri)cos  Qg~R;i(p-R-0 
Rl  R3  sin  (J  2 

R3  lo231ii39UO 

(P-fil)  ^  -Oo3862lil39 
cos     02  0.9915U671 

R-,  1.216791*93 

(p-R3)  -0.1*0088586 

R3(p-R1)  .  -O.U7563287 

R3(p-Rl)cos    02  -0.1*7161221 

R1(p-R3)  -O.U8779588 
Rjup-Rijcos  Cf2~Rl(P™R3)  O.OI6I8367 
sin  (J  2  0.12975025 

R1R3  sin    (J '  0.191*1*1897 

e  sin  (ui  -  4)  0.08321*121 

e  cos   (ui  -  w)     a  (p-Ri)/Rl     -     -0.3171*2521 

c2sin2(Ul«w)  0.00692910 

e2cos2(ui-w)  0.10075876 

e2[sin^(ui-w>fcQs2(ui-w)l  0»  107657  86 
e   L  J      0.32815828 

p/l-e2     .  a  0.93078812 

e  sin  (ui=w)  0.08321*121 

sin  (ui-w)  0.25366177 

e  cos  (ui-w)  -0.3171*2520 

cos  (ui-w)  -0.96729298 

(ui-w)  165°18«20".508 

ui  2U8013502».8U5 
w  82o5U»li2».337 


V  *  u-w 
V/2 


1-e 
1+e 

V(l-e)/(l  +  e) 
tan  V/2 
tan  E/2 
E/2 


E  (radians) 
sin  E 
e  sin  E 
M  (radians) 
M 


TABIE     10 
COMPUTATION  OF  TIME  OF  EPOCH,  T0 
Observation  #1  Observation  #2 


lil 


165°18»20M.508 
83o39«10».25ii 


170°19907".97U 
85°09B33,t.987 


tan  E/2  =   "Y("l«e)/(1+ e)  tan  v/2 


7.7555753 
5.5159787U 
79°U3'27".723 
l59°26»55l,oUli6 


0o67l8Ul72 
1.32815828 
0.5058UU62 

0.7U2275U 

llo808U28 
8o398U7920 
83°12«35"«307 

l66o25810».6lU 


M 


E  -  e  sin  E 


2.782905U- 
0.3510U533 
0.115198U3 
2.66770698 

153°25816».7H 


2.90U57022 
0.23^80935 
0, 077051*63 
2.82751559 
162000 « 16" o 966 


n    =    k/VaT 


k  *- 

a3  *" 

v^ 

n  9- 

E-e  sin  E/nsb^T0  0.02237350 
tj  (days) 

j  =  (1,2,3)  Dec  18.93692292 


107.0731376 
O0806U0366 
Q. 89799981 


Observation  #3 

172°ii5839"o087 
86°22«U9".5i43 


I5<.80837it 
U.2U335095 
8U°55»02».65U 
169°50«05!,.308 


2.96U7658 
O.I76U8680 
0.05791560 
2 . 90626098 
I6603C  59  ".358 


119.23511*50  r  ado /day -* 

0.02371378  0.02U37li20 

Dec  I8.9382li560        Dec  18.93892U768 


To 


Dec  18.91U5U9U2        Dec  18.911*53182        Dec  I8.91ii55056 


TABLE     11  U2 

ORBITAL  ELEMENTS  DETERMINED  FROM  PUBLISHED  DATA 


(t-T0)  *    -U.085U6818  days  w  •  110°  U8»l*0«o896 

sin  w=    0.93U7&25 
(t-T0)2*     16.6910^025  cos  ws-0o35529231 

(t-T0)3-  -68.19075U69 

w      s  (li+0°a6)  +  (7°.275)(»ii«>0851i68l8)-|-(0.OOOU27)(l6a69105025) 

-(°.22U)(-.3#29231) 

w      =  110°  $V   $3"*l$h 

Q     =  (l86o.278)~(Uo.896)(«lio08^68l8)-(o.000298)(l6<)69105025) 

-(O.00887)(-«3#29281) 

Q     =  206°  16«  li3".067 

i   =  (330.210H0oOOO^U)(-Uo08^68l8)  +  (°o000005)(l6o6910$025) 

+  (0oOOii73)(o93U75^) 

i   -  33°  13 f  00».U02 

e   =  (.0997lH.000202)(«.U<,08^U68l8)+(o00000003)(l6o69105025) 

~(.OOOU^19)(o93li7#2£) 

e   =  0.10011028 

a   =  q/l-e  =  1.0^3706/„89988972 


U3 

Comparing  the  above  results  it  is  apparent  that  the  position  of  the 
orbital  plane  is  correct  but  that  there  is  a  discrepancy  in  the  orienta- 
tion of  the  orbit  in  its  plane  and  in  the  size  and  shape  of  the  ellipti- 
cal orbit o     Reviewing  the  computations,  it  can  be  seen  that  the  auxiliary 
parameter  p,  which  has  a  value  of  approximately  0,830,  is  too  small.     In 
the  case  of  an  elliptical  orbit  it  should  be  slightly  greater  than  one0 
Since  p  is  a  function  of  the  auxiliary  parameters  uj,  it  is  believed  that 
the  source  of  the  discrepancy  rests  here0     This  would  affect  the  results 
in  the  computations  for   (j*  and  for  the  argument  of  the  perigee,  w, 
also.     To  avoid  continuing  the  problem  with  this  discrepancy  present  I 
decided  to  use  the  published  orbital  data  from  the  82  observations  rather 
than  to  continue  with  the  non-realistic  values  for  e  and  a  derived  from 
the  computations.     The  orbital  parameters  computed  for  the  times  of 
observations  in  the  next  chapter  are  based  on  those  published  in  Smith- 
sonian Astrophysical  Observatory  Special  Report  #78 • 


£•  DETERMINATION  OF  THE  GEOGRAPHIC  COORDINATES 
OF  THE  UNKNOT/N  STATION, 


There  are  several  methods  available  to  determine  the  position  of  an 
unknown  station  from  observations  of  a  celestial  body.  Two  of  these 
methods  were  devised  by  William  Markowitz  of  the  U.  S.  Naval  Observatory 
for  his  well  known  Moon  Camera  Method  and  can  be  adopted  for  use  with 
artificial  earth  satellites.  These  methods  are  called  the  linear, 
which  gives  the  coordinates  of  an  unknown  station,  and  the  differential, 
which  provides  a  way  to  find  the  corrections  to  an  assumed  approximate 
position.  The  differential  method  has  been  selected  for  use  in  this 
particular  problem.  As  mentioned  earlier  in  this  paper,  only  two 
approximations  in  the  differential  method  will  be  made. 

The  procedure  to  determine  the  coordinates  of  the  unknown  station 
is  as  follows: 

1.  The  orbital  elements  of  the  satellite  are  known  for  a  particular 

time,  T0,  the  time  of  epoch.  Now  assuming  these  elements  are  correct, 

-  .   i 

the  orbital  elements  are  then  determined  for  the  times  of  the  observa- 
tions from  the  unknown  station. 

2.  Using  these  elements  the  geocentric  right  ascensions  geocen- 
tric declination,  and  the  range  of  the  satellite  from  the  earth's  center, 
will  be  predicted  for  the  times  of  the  observations  from  the  unknown 
station. 

3.  The  differences  in  the  computed  coordinates  and  the  observed 
coordinates  of  the  satellite  will  furnish  us  with  a  set  of  coefficients 

for  three  differential  equations  to  be  solved  for  each  observation.  The 

i 

solution  of  the  differential  equations  will  produce  the  increments  IXCC 

c  ' 
and  A  0       which  are  added  to  the  approximate  coordinates.     Using  these 


new  coordinates,  the  differential  process  will  be  repeated  and  the  resul- 
tant increments  will  be  added  to  the  approximate  coordinates  again.  This 
will  be  the  final  result. 

Refer  to  Figure  6, 

The  observer,  0,  views  the  satellite  along  the  line  OS'   at  a  distance 
of  r.     The  coordinates  observed  from  this  position  are  the  apparent  right 
ascension,    (X   ,   and  the  apparent  declination,    0     .     From  the  center  of 
the  earth,  E,  the  satellite  would  be  seen  along  the  line  ES'   at  a  range 
of  R.     The  coordinates  observed  from  this  position  are  the  geocentric 
right  ascension,  CX  s  and  the  geocentric  declination,     0    o     The  coordi- 
nates of  the  satellite  derived  from  the  orbital  elements  are  the  geo- 
centric coordinates  and  will  differ  from  the  observed  coordinates  determined 
from  the  camera  film  by  the  amount  of  arc,  SS',  in  the  orbital  plane. 
This  difference  enables  us  to  determine  the  observer's  position  on  the 
surface  of  the  earth.     This  difference  between  the  computed  and  the 
observed  positions  of  the  satellite  should  be  very  small  but  not  so 
small  that  the  two  positions  would  be  at  the  zenith  of  the  observer .     At 
this  position  both  the  computed  and  observed  positions  would  have  the 
same  value  and  the  observer's  position  could  not  be  determinedo     Another 
indeterminate  situation  would  exist  if  the  satellite  was  observed  at 
the  horizon.     Here  the  parallactic  angle  would  be  too  large  and  the 
effect  of  refraction  would  make  the  observations  uselesso      [l2,  p.3lu  • 

Proceeding  with  the  problem,  the  following  data  are  knowns 

Location  of  Observer:      j.  Curacao,  N.W.I. 
Approximate  Positions     ?\     -  12°  13 ' 23" • 21  N 

X    -  68°  55* 07". 86  W 

xa  -      0.35199281*2 

ya  *       0.913105082 

Za  «         0.210569H7 

Satellite  Observed %  Explorer  I  (1958  Alpha) 


U6 


hi 

Universal  Time  of  Observations:     23h5lrak7s.OU 

23&1W.98 
23h$2m12Sol£ 

The  approximate  rectangular  space  coordinates,  xa,  ya,  za,   are  found  in 

the  same  manner  as  that  used  to  find  the  coordinates  of  the  known  station, 

This  procedure  is  discussed  in  Chapter  k  and  will  not  be  shown  here. 

Using  the  above  data,  the  orbital  parameters  for  the  times  of  obser- 
vations are  computed.     These  computations  are  contained  in  Table  12„ 

With  the  orbital  parameters  known,  the  next  step  is  to  compute 
the  eccentric  and  true  anomalies  of  the  satellite  using  this  formula 
[U,p.l6lj    for  the  eccentric  anomaly: 


E  ■  M+e  sin  M.    ®£  sin  2  M  + 


o  O        O   o    o 


where  M  is  in  radians  <>     The  additional  terms  in  this  formula  are  not 
listed  or  used  in  the  determination  of  E,     Moulton  states  that  because 
of  the  small  eccentricity  of  the  orbit  the  series  converges  very 
rapidly  and  the  first  three  terms  give  an  accurate  value  of  EJo 
Then  the  true  anomaly  V  is  computed  from  the  formula: 


tan  |  V  =  ~V(l-f-e)/(l-e)  tan  | 

The  computations  of  the  anomalies  for  the  three  observations  are  listed 
in  Table  13 „ 

The  next  step  is  to  determine  the  geocentric  right  ascension,  geo- 
centric declination  ^  and  geocentric  distance  of  the  satellite  fbr  each 
observation  using  the  formulae s 

sin  U     -  sin  i  sin  (w+v) 
OC  *    Q  +    arctan    [cos  i  tan  (w  +  v)J 
R    s  a  (1  -  o  cos  E) 
where  a  is  the  semi-major  axis  of  the  elliptical  orbit  in  terms  of  the 


TABLE  12  U8 

OBSERVATION  1   OBSERVATION  2   OBSERVATION  3 


(U.To) 


Location 

Time  of  Observation 

Date  of  Observation 

Smithsonian  Date  of 

Smithsonian  Date  of  Epoch  (Tn) 

(t-T0) 

w  (Mean) 

sin  w 

cos  w 

(t-To)f 

(t-T0)3 


Curacao,N„W.i 
lffil*7s.0l* 


23% 

18  Daeanber  1560 

Observation^ )  37286.991*291*1*1*1* 

37291.0 

-U.005705556 

111°23,27".7255 

.93U1290 

-.361*73110 

16.01*567700 

-6Uo27U2575l 


CuracaOjN.W.i. 
23h5lm51*s.98 
18  Bsaanber  I960 
37286. 991*38631*3 

37291.0 

-U.005613557 

Hl°23'30».1299 

.93110865 

-.361*71*195 

16.01*1*91*077 

-61*.  26983387 


Curacao  ,N  .W.I. 

23h52m12s.l*50 
18  DBcemb3rl960 
37286.991*58651*2 

37291.0 

-I*.o651*ll!*5o 

lllo23«35,,Ji2  00 
,93109929 

-.361*76583 

16.01*332095 

-61*. 260101^6 


»  (ll*0°.6   ±2)+70.275(t-To)+00.OOOl*27(t-To)2+00.22l* 


cos  w 


7°.275(t-T0) 
0°.  0001+27  (t-T0)  2 
0°.22l*  cos  w 

wto 

*t 

»t 


-29°oll*l50792 

0°.00685l50 

-0°. 08169977 

12*0°. 6  ±2 

Hl0.38361*38l 

llio  23«01M.ll8 


~29°.U*083935 

0°. 00685119 

-00.08170220 

U*0°.6  ±  2 
1110. 381*30961* 
lllo23'03",^ 


-29°.13936836 

0°. 00685050 

-0O.08170755 

U*0°.6  ±  2 
lll°.38577ii59 
Hlo23'08Mo789 


Q=  (I860.278  ±9)  -  (itO.896(t-To)-0°.000298(t-To)2+0o.00887  cos  w 


U°.896(t-T0) 
0°. 000298 (t-T0) 2 

0°*DO887  cos  w 


19°.6ll93l*l*Q 
-OO.OOU78161 
-00.00323516 
1860.278  ±9 
205O.88191763 
205°52  «51iM<?03 


19°.6lll*81*l*6 
-0°. 001*78139 
-0°.00323526 
186°.278  ±9 
205°.  8811*6781 
205o52'53H^81* 


19°.6l0U9l*50 
-0°.  001*78091 
-0O.003235U7 
186°. 278  ±9 
205°.  8801*7812 
205°52'!*9IV721 


i  =  (33°.210  ±3)  ^o.00055i*(t^o)  +  ^oo^CxlC-*^(t-To)2-<jO.00l*73  sin  w 


0°.000551*(t-To) 
O°.000005(t-T0)2 

0°.OOl*73  sinw 
i  To 

it 


0O.00221916 
0°o00008023 
-0°.00l*l*0l*l6 
33°«»210  ±  3 
33O.20789523 
33°12»28'»Jj25 


00.00221911 
0°o00008022 
-0°.00l*l*0l*ll* 
33°.210  ±  3 
330.20789519 
33012'28n.l£3 


00.00221900 
0°. 00008022 
-0°.  001*1*01*10 
33°.210  ±3 
33° .20789512 
33°12'28M.1*22 


e  s  (.09971±8)-Oo000202(t-To)+Oo30xlO-7(t-T0)2+ 0.0001*5519  sinw 


0. 000202 (t-T0) 
0.00000003 (t-T0)2 

0.0001*5519  sin  w 

©To 

«t 


0o000809l5 
0.0000001*8 

0. 0001*2383 

0.09971 ±8 
0.10091*31*6 


0.00080913 
0.0000001*8 
0.0001*2383 
0.09971+8 
0.10091*31*1* 


0.00080909 
0.0000001*8 
0.0001*2383 
0.099711  8 
0.10091*31*0 


(continued  on  page  1+9) 


TABLE  12  (Continued) 
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OBSERVATION  1 

OBSERVATION  2 

OBSERVATION  3 

M«  (.9li0l6)  +  (13.]*35l3U)(t-ToH0.0002839(t-To)2  -0o  11x10"^ (t-T0)3- 

-0.000631jc10~3 
cos  w 

(13.l*35l3U)(t-T0)o 

-53 0 81719091 

-53.81595623 

-53.8I323966 

(0.0002839)(t-To)^ 

O.OOU55537 

O.OOU55516 

O.OOU55^70 

(0.0000011) (t-T0)3 

0.00007070 

0.00007070 

O.OOOO7069 

(O.OOOOOO63)  cos  w 

0.00000023 

0.00000023 

0.00000023 

MTo 

0.91*016 

0.9U016 

0.9U016 

M 

-53«8l256U6l 

-53«8ll3302it 

-53«8o86lliOU 

ttb  (Revolutions) 

0.12759539 

0.12882986- 

0.1315U596 

*t 

U5°.93U3UOUO 

U6°.3787li9(SO 

ii7°c3565U560 

*fc 

U5°56'03"j625 

U6°22>li3n<Jj99 

U7°21«23".5& 

TABLE  13 
COMPUTATION  OF  ECCENTRIC  AND  TRUE  ANOMALIES 


50 


OBSERVATION  1 

OBSERVATION  2 

OBSERVATION  3 

S  =  M  +  e 

e2 

sin  M+*2*     sin  2  M 

e2 

0.01018958 

0.01018958 

0.01018957 

e2/2 

0.005091*79 

0.00509U79 

0.005091*79 

1  +  e 

1.10091*31*6 

1.10091*314* 

1.10091*31*0 

1  -  e 

0.89905651* 

0.89905656 

0.89905660 

(l+e)/(l-e)  1* 

1      sin  M        J 

1.22U55I*19 

1.221*551*11* 

1.221*551*01* 

1.10659577 

1.10659571* 

1.10659570 

0.71851*327 

0.72391601* 

0.73558351 

sin   2M 

0.9991*6819 

0.99881*211 

0.99661861* 

e  sin  M 

0.07253221* 

0.073071*58 

0.071*25230 

(e2/2)sin  2  M 

Oo 00509208 

0.00508889 

0.00507756 

E  (radians) 

0.87932980 

0.88762535 

0.90585750 

E 

50?3818863U 

50985718631* 

51990181159 

E 

5O°22'51+H.7S0. 

50°5lt25,,.8?l 

5l°51*I06».522 

E/2 

25°11,27,,.396 

25°25'1*2».986 

25°57»03».26l 

„tanE/2        - 

0.1*7037122 

0.1*751*1*692 

0.1*8667231* 

|{1^« )/{!-•> WtmS/ZtstR  V/2 

L        v/2         J 

0.52051080 

0.5261275U 

0.53851*952 

27029'50».872 

27°1*5'00».338 

28018' 16».789 

V 

51*059'l*l».7l*l* 

55°30«00».676 

56036' 33".  578 

TABLE  Ik 
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COMPUTATION  OF  GEOCENTRIC  DECLINATION,   GEOCENTRIC  RIGHT 
ASCENSION,  AND  GEOCENTRIC  DISTANCE  OF  THE  SATELLITE-, 

OBSERVATION  1       OBSERVATION  2       OBSERVATION  3 

sinO  =  sin  i  sin  (w+V  ) 


(w+v) 
sin  (w+V  ) 
sin  i 
Bin  (5 

6 


166°22»U2»o862 
0.23550558 
0.5U767853 
0.12898135 
7°2U838».636 


I66053*0li».191 
0.226911*83 
0.51*767853 
0.121*27638 
7O08 « 20". 290 


a  = 


tan  (w+v  ) 

cos  i 

cos  i  tan  (w+v  ) 

Q       r  ,      0 

arctan    cos  1  tan(w+v  ) 

a     L 


q  (Ha) 

q  (earth's  radius) 


W+  arctan  [cos  i  tan  (w+v  ) 

-0.2U232136  -0.2329925U 

0.83668886 
-0.20271*758 
205°5285U"«903 
168032" 19" .601 

U^'lU^Ol* 


O.83668886 
-0.191*91*226 
205°52«53"o281i 
168°58»08».322 

U*o5l»01".6o6 


E 

cos  E 

e  cos  E 

1  -  e  cos  E 

R 


a  =  q/  (1-e) 

6.72127351 
1.05375739 
1.17207021 

R  ■  a  (1  -  e  cos  E) 

0.10091*31*6 
50°22«5li".791 
0.637667& 
0.061*36837 
0.93563163 
1.09662596 


6.721273^9 
lo05375739 
1.17207018 


0.10091*31*1* 
50o5l«25".871 
0.63125552 
0.06372110 
0.93627890 
1.097381*58 


167°59,U2".367 
0.20799531 
0.5U767853 
0.113911*57 
6°32«27".682 


-0.21261*591 
O.83668886 
-0.1779181*6 
205°52'U9".721 
l69°51i,Ul"o720 
l50U7»3l"Jil 


6.7212731*9 
1.05375738 
lol7207012 


0.10091*31*0 
5l°5U806».522 
0.61701099 
0.06228319 
0.93771681 
1.09906985 


*2 

earth's  radius.     The  value  "a"  can  be  computed  from  the  formula: 

a  =  q/(l-e) 
where  q  is  the  distance  of  the  satellite  at  perigee.     The  computation 
of  these  values  are  shown  in  Table  lit.      C/g,  the  Greenwich  hour  angle 
of  the  vernal  equinox^,   is  tabulated  in  Table  l£0     The  computations  of 
the  sidereal  times  of  observations  from  which  this  value  is  derived  is 
not  shown  or  listed.     The  sidereal  time  was  established  in  the  same 
manner  as  it  was  for  the  observations  from  the  known  station. 

Now  the  computed  apparent  right  ascension  and  declination  must 

be  found  using  the  formulae  [12,  p.37J    : 

x        t  /-v       rv  \  ,      xa  sin  G  +  y,  cos  G 

tan  (CXc-  (X  )  ae_4cosg^y^sinG.fftcos($ 

and 

tan      6«     =      cos  (CCqrCC)  (R  sin  6     -a) 

7      "  -Lxa  cos  G  +  ya  sin  G  +R  cos  X 

Then  the  differences  between  the  computed  and  the  observed  values 

are  determined  byr 

A  a'  =  a'-  ocl 
A  6'  =  6-  6c 

The  values  and  computations  of     0  c  and     OC  c  and  the  increments 
A  OC     and   A   0      ape  tabulated  in  Tables  IS  and  16. 

The  differential  equations  to  be  solved    [l2sp.36j     ares 

i 
ai  dxa  +  a2  dya  s  a$  d  QC 

b!  dxa  +  b2  dya+b3  dza  '=  b^  dCX  +  b£  dO 

The  coefficients  of  these  equations  are: 

ax  ■  sin  G  +  tan  (  C£  -  OC )  cos  G 

a2  -  cos  G  -  tan  ( (X  -  CO  sin  G 

a3  =  sec2  (CX-CX  )(R  cos  6    -  *a  cos  G  +  ya  sin  G) 


53 
b^  =  tan  0    cos  G 
b2  =  -  tan  0     sin  G- 

b3  =  -  cos  ( a'  -  a  ) 

b|j  s  sin  (  (X-  (X  )     (R  sin  0  -  za) 
bs;  ■  sec2   0   (R  cos  0  -xa  cos  G  +  ya  sin  G) 
The  computation  of  these  coefficients  are  contained  in  Table  16.     The 
solution  of  the   first  approximation  differential  equations  is  shown  in 
Table  17.       ACIC  Technical  Report  No.  86  suggests  these  equations  be 
solved  by  Doolittle's  method  but  being  more  familiar  with  the  Gaussian 
solution  of  normal  equations ,  I  decided  to  use  the  latter  method*     The 
solution  of  the  normal  and  reduced  equations  yields  the  increments 

A  xa».    A  ya*  and  Aza,     These  increments  are  then  added  to  the  origin- 
al approximate  coordinates  to  find  the  corrected  position  of  the  observer, 
The  corrected  space  coordinates  ares 

x0  y0  z0 

+  0.3£LU627£3  +  0.912606236  +  0.21092082U 

Since  the  International  Ellipsoid  was  used  as  the  reference  ellip- 
soid,the  xa,  ya,  za  coordinates  can  be  multiplied  by  6,378,388  meters, 
the  semi-major  axis  of  the  ellipsoid,  to  find  the  metric  distance  of  r , 
the  radius  of  the  observer's  position  from  the  earth's  center.     This  has 
not  been  done  lecause  it  is  of  little  interest  in  this  problem.     How- 
ever, the  distance  of  the  observer  from  the  earth's  center  is  found  in 
units  of  the  earth's  radius  by  the  equations 

P  s_v/xo  +  To  +z0'       »  1.000U31808 
Then  the  resultant  geographic  coordinates  of  the  observer  may  be  easily 
derived  as  follows: 

sin  (jf/  =  z°/p     -  +  0.21082979 
tan    0'       =  +0.21^67762 


TABLE  15 

/(APPARENT  RIGHT  ASCENSION,  QCc I 
SOLUTION  FOR  COMPUTED  V APPARENT  DECLINATION,  ff 


FIRST  APPROXIMATION 


51* 


CC 

6 

R 


8 


9g  (gha/) 

oc-  9g*Q 

sin    G 
cos  G 
sine 
cos  ^ 

(  ©:'-  ex) 


sin  (CC 
cos  (  GC 
tan  (^ 
sec2(C(' 
tan.  6' 
see2  £' 


-CC) 

-a) 
-co 


OBSERVATION  1 

7°2li,38».636 
1.09662596 
355°  1*6«30". 
-30°  5l800««. 

&5°31°51*,,.975 

-71°06»  1*0". 1*71 
-0„9l46ll;889 
0.32373178/ 
0.12898135 
0.99161*703 
3i*l°21«l5M.l*96 
-0.31971509 
0.91*751372 
-0.3371*2529 
1.11385581* 
-0.59730303 
1.35677091 


tan  (CXC  -a) 


OBSERVATION  2 

lU°5lt01».6o6 
7°08'20,,.290 
1.097381*58 
359°  39'  CO". 
-32°  52 •  00". 
85°33«51*n.l*o5 

=70°l*2«52tt.799 
-0.91*388553 

0.33027279 

0.121*27638 

0.99221*765 
3l*i*°l*7,'58«.39l* 
-0.26219670 

0.96501U95 
-0.27170236 

1.07382217 
-0.61*6101*06 

1.1*171*501*6 

xa  sin  G-fya  cos  G 

=oca  cos  G  -J-  ya  sin  G  +  R  cos  Q 


OBSERVATION  3 

l5°li7'31"oUil 
6032* 27" 0682 
1.09906985 

8D  03'  00". 
-36°  39'  00". 
85°38»l7«.l60 

-69°50'1*5"»719 

-0.93877011* 

0.31*1*51*1*07 

0.113911*57 

0.9931*9055 

-7°UU'3r'oUa 

-0.131*71371* 

0.990881*56 

-0.13595301 

1.0181*8322 

-0.71*1*0201*0 

1.55356636 


Aa 

ya 

za 

xa  sin  G 

ya  c°s  G 
-xa  cos  G 

ya  sinrG 
R  cosQ 


xasinG+ya  cos  G 
-xacosG+yasinG+ 
4  cos  () 


tan  ((Xc-CO 

(a'c-  a) 
a' 


-0.3330376U 

0.29560113 

-0.11395127 

-0.86393336 

1.0871*6588 

-0.037143651 
0.10958125 


-0.31*16321*1* 

3l*l008«17M.l*06 

355°33»31wo910 


tan 


6', 


cos 


+  0.35199281*2 
+  0.913105082 
+  0.210569117 

-0.33221*095 

0.30157376 

-O.II625366 

-O.86I86667 

1.08887727 

-0.03066719 
0.11075691* 


-0.27688730 

3l*l*°31«23M.757 

359°22«25".363 


-0.3301*1*037 

0.311*601*91* 

-0.12127705 

-0.85719579 

1.09191551 

-0.0158351*3 
0.1131*1*267 


-0.13958972 
-7°56«1*7».593 
70501  itf"  .81*8 


((*'    -00  (R  sin6-za) 


°xa  cos  G+ya  sin  G+R  cos 


TABLE  15  (Continued) 
FIRST  APPROXIMATION 


& 


OBSERVATION  1 

cos  (CLq-OL)  0.9U630093 

R  sin  Q  Q.1)|1)|)ilt30 

(R  sin  0  -za)  -0.06912U82 

-xa  cos  G  +  ya  sin  G  + 

R  cos  fa  0.109^8125 
oosCa'    f-CO* 

(R  sin  0-za)  -O.O65I4I288 


OBSERVATION  2 

0796373889 

0.13637898 

-0.07U190114 

0.11075691* 

-O.O7H49992 


OBSERVATION  3 

0.990397U6 

0.12520007 

-0.08536905 

0.113Uii267 

-0,081i5it929 


tan 

6! 


-0.59693U97 
-30°50'OU".037 


-0.6U555702 
■32°50'iiO,,.376 


-0.7U530U13 
-36°U1'50»,335 


TABLE  16 
COMPUTATION  OF  COEFFICIENTS  OF  DIFFERENTIAL  EQUATIONS 
FIRST  APPROXIMATION 
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OBSERVATION  1 


OBSERVATION  2 


&i  ■  sin  G  +  tan  (  (X  -  CO  cos  G 


a-ac  = 


6 


•& 


(X 
(radians) 

A5A6< 


0°13'58».O9O 
0°. 232 80278 
O0OOU063I8 
0°01«OUM.037 
0°. 01778806 


^^'(radians)  -O0OOO3IOU6 


sin  G 

tan((X-O0cos  G 
al 


cos  G 


) 


tanCa^OOsin  G 
*2 


sec2((X'-CO 
R  cosQ 
xa  cos  G 
ya  sin  G 

RcosQ-xa  cos  G+ 
ya  sin  G 

»3 


sin  (ff'-OQ 

R  sin  Q 

za   /: 

R  sin  0  -z_ 


-Oo9U6LU889 
-0.10923529 
-1.05538U8 


0°l6'3liM.637 
0° 027628805 
OoOOU822lU 
-0°01«19".62l, 
-0°. 02211778 
-0.00038603 

-0.9U388553 
-0.08973590 

-I.03362H43 


OBSERVATION  3 


0°12'l6".l52 

0°.20itU8666 

0.00356897 

0002 «^0» .335 

0°.0U731527 

0.00082581 

-0,93877011; 
-O.OU68U180 
-0.9856119U 


a2    =  cos  G  -  tan  (QC-  GO  sin  G 


0.32373178 
0.31925U56 
O.OOW7722 


0.33027279 
0.256U5593 
0.07 3 81686 


0.3l^5UU07 
O.12762863 
0.2l69l5Wi 


a$  =  sec2  (CX  -O0(R  cos(5-xa  cos  G+ya  sin  G 


1.1138558ii 

1.087U6588 

0.11395127 

-0.86393336 

0.10958125 
0.12205772 


lo07382217 

1.08887727 

0.11625366 

-O.86I86667 

0.1107569U 
0.11893326 


b]_    —  tan  0       cos  G 
-0.19336597  -0.21339059 

D2    =  -  tan  U  sin  G 
-0.56513760  -0.6098U827 

b3    =    -  cos  (CX'-OO 
-0.9^751372  -O.965O0M5 

b^  =  sin  (Ot'-CX  )(R  sin  (5  -za) 


-0.31971509 

Q.lUl)iMi30 
0.21056912 

-0.06912U82 

0.02210025 


»0o26219670 

0.13637898 

0.21056912 

•0.07Ul9011i 

0.019U52U1 


1.018U8322 
1.09191551 
0.12127705 

-0.85717579 

0.113UU267 
0.115539U6 


-0.2563U782 


-0.698U6Uli* 


-0.99088U56 


-0.13U7137U 
0.12520007 
0.21056912 

-0.08536905 
0.01150038 


sec 
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TABLE     16     (Continued) 

FIRST     APPROXIMATION 
OBSERVATION  1  OBSERVATION  2  OBSERVATION  3 


b£    =  sec2  6    (R  cos  Q  -  xa  cos  G+ya  sin  G) 

,w2(§                         1.35677091              1.10.7U50U6  1.55356636 
R  cos  Q    -  xa  cos  G+ 

ya  sin  G            0.10958125                0.11075691*  0.1131^267 

b£                                   0.11*867  665                0.156992U8  0.17621072 


TABLE  17 
SOLUTION  OF  DIFFERENTIAL  EQUATIONS 
FIRST  APPROXIMATION 
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«,  Aa 

b*A6' 
bUA«/+b*A6' 


OBSERVATION  1 

0.0001*9591* 

0.00008980 

-O.OOOOl*6l6 

0.00001*361* 


OBSERVATION  2 

0.00057351 

0.00009380 

-O.OOOO606O 

0.00003320 


OBSERVATION  3 

O.OOOU1236 
0.00001*101* 

0.00011*551* 
0.00018658 


dxa 

-1.055381*18 
-0.19336^97 
-1.0336211*3 
-0.21339059 
-0.98561191* 
-0.25631*782 


dxa 
3.302280071 


OBSERVATION  EQUATIONS 


aidxa  +  a2<iya  =  33  &CC1  .  r 

b^a  +  b2dya  +  b3  dza  =  b^  da'    +    b£d6 

d7a 


ciz. 


0.001*1*7722 

-0.56513760 

0.07381686 

•0.60981*827 

0.2169151*1* 

=0.6981*61*11* 


dya 

0.12361*5699 
1.231668857 


-0.91*751372 
-0.965011*1*5 
-0.990881*56 


-0.0001*9591* 
-0. 00001*361* 

-0.00057351 
-O.OOOO3320 
-0.0001*1236 
-0.00018658 


NORMAL  EQUATIONS 

dza 

0.61*3153009 
1.816085355 
2.810887350 


0.001585979 
0.00001*1227 
0.000258267 
0.000000783 


-1.0511*0290 
-1.70606093 
-0.96037808 
-1.78828651 
-0.76910886 
-1.91*588310 


1*„070661*758 

3.171MI1137 
5. 2703 83 981 
0.001886255 


^a 

dxa    3.302280071 

dya  -0.0371*1*2523 

dza  -0;i9l*760285 

-0.0001*80268 


dya 

0.12361*5699 

1.227039250 

-1.1*601*29309 

0.000011*797 


REDUCED  EQUATIONS 

dZQ 


0.61*3153009 
1.792001*081* 
0.0685311*01 
0.000351707 


dxa  dya  dza 

-0.000530089     -0.0001*9881*6  +0.000351707 


0.001585979 
-0.000018156 
-0.000021*103 

0.000000013 


1*.  070661*758 
3.019025178 
0.068507298 
0.000000012 


xa        0.35199281*2 
Axa      -P.  000530089 

xo  +  0.3511*62753 


ya       0.913105082 
l\ya    -0.0001*9881*6 

y0+  0.912606236 


za        0.210569117 
Aza         Oo  000351707 

z0+  0.210920821* 
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tan  0=  tan  0/l-e2     =    0.21567762/0.99327733    =  0.21713736 
tan  A      =      7o/Xo    =  2  0^965910*6 
which  result  ins 

0    =    12°  15 «  03"  ,168  N 
A     =    68°  56'  Hitto336  W 
The  differential  method  was  repeated  again  for  a  second  approxima- 
tion.    The  corrected  space  coordinates  from  the  first  differential 
process  were  used  as  the  approximate  observer's  position.     The  data 
used  in  and  the  computations  of  the  second  differential  approximation 
are  listed  in  Tables  18 ,  19,   and  20. 

The  corrected  and  final  results  of  the  second  approximation  are: 

x0  y0  z0 

+  0.35Hi8350U  +  0.912709058  +  0.21081*9193 


P    =  ~\j*l  +  7o2  4  zo2'     =  Io0005l7795 
sin  0';  =  zo/p     =     0.2107U007 
tan  0'  =     0.21558157 
tan  0  =   tan  0/l-e2  =  0.21558157/0.99327733    =  0.21701*066 
tan  A   =  yb/xo  =  2.59673370 
The  corrected  and  final  geographic  coordinates  are: 
0    =    12°  11a«  W.118  N 
A    =    68°  56 »  18».0U5  w" 
Further  approximations  would  probably  yield  changes  to  these 
coordinates  but  as  stated  previously  the  assumption  that  two  approxi- 
mations are  sufficient  has  been  made. 

As  a  note  of  interest^  and  to  satisfy  my  own  curiosity,  further 
calculations  were  made  using  the  linear  method  but  are  not  included  in 
this  paper.     First,  I  used  my  own  computed  value  of  the  Gaussian  gravi- 
tational constant,  k,  in  the  formula  to  determine  "a": 


TABLE    18  „  60 

S0UJTI«  op  c«h,m  f^sTT  ss:™-  gj 

SECOND  APPROXIMATION 
OBSERVATION  1       OBSERVATION  2       OBSERVATION  3 

tan  (  OCc-CX  )    _  *a  sin  G  +  ya  cos  G 

-xa  cos  G  +  ya  sin  G  4-  R  cos  £ 


Aa 

7a 
za 

xasin  G 
yacos  G 
-xa  cos  G 
ya  sin  G 
R  cos  fa 

xasinG+yacosG 
-xacosG-tya  sin  G  + 
R  cos  (5 


tan  (CXc-O0 
(OGrOO 


cxc 


-0.33253609 

0.295U396U 

-O.U377966 

-O.863U6138 

1.087U6588 

-0.037096U5 
Ooll022U8U 

-0.336^251i 
3ia°23857"ol23 
355°U98ll"o627 


+  0„3^1i462753 
+  0.912606236 
+  Oo21092082U 

-0.3317U060 

0.3011*0901 

-0. 11607858 

-0.86139582 

1,08887727 

-0.03033159 

0.1111*0287 

-0.27226938 

3iiit°U6«09tt.U92 

359037'11"<098 


tan  A'     _     cos(a'c-CX  )(R  sin  6  -za) 


xa  cos  G4-ya  sin  G+R  cos 


cos  (OCc-OQ 

R  sin  Q 

R  sin  8    -za  r 

cos(Q^-00(R  sintUa) 

•-54  6<c 


0.91*776396 
0.11*11*1*1*30 
-0.0691*7652 
-0.06581*731* 
-0.59739111 
"30°5l813"o391 


0.961*87589 
Ool3637898 
-0.071*51*181* 
-0.07192362 
-O.6U561730 
-32°50'l*9M.l52 


-0.32991*271* 

0.311*1*3307 

-0.121091*1*1 

-0.8567271*9 

1.09191551 

-0.01550967 

0,llli0936l 

-Ool35938ll 
-7°l*i*828M.i*2l* 
8°  03 8  03  ".017 


0.99088653 
Ool2520007 
-0.08572075 
-0.081*93951* 
-0.71*1*1*7237 
-36°39859".995 


TABLE  19  61 

COMPUTATION  OF  COEFFICIENTS  OF  DIFFERENTIAL  EQUATIONS 
SECOND  APPROXIMATION 
OBSERVATION  1   OBSERVATION  2   OBSERVATION  3 


ax  =  sin  G  +  tan  (CC-CQ   cos  G 


A  OC  (radians) 

6'-6:=A6; 
„/     A<5' 

^Q  (radians) 
al 


a3 

sec2  (OC-CO 
R  cos  Q 
xa  cos  G 
y.  sin  G 

(*  cos  (5-3ca  cos  G  + 
ya  sin  G) 


sin  (QC'-OO 
R  sinQ 

Z*   X 

R  sin  0  -  za 


-0°02»UlMo627 

-0°oOW+81i639 

-0.00078272 

0°00«13»o391 

0°o00371972 

Oo000061i92 


0°01Bli8»,912 
0°. 0302^333 
0.00052802 

-0001»10,'<>8ii8 
-0° .01968000 
-0. 0003143  81 


a2  = 


-lo05538Ul8    -1.0336211*3 
cos  G  -  tan  (Ot'-OQ   sin  G 
0<,OWi7722     Oo0738l686 


-0°00«O3tt«O17 

-0°o00083806 

-OoOOOOli63 

0°00» 59 "o995 
0°. 01666^28 
0.00029086 

-0.9856119U 


0o2l69l5U; 


=  sec^  (CX  -O0(R  cos  6  -  xa  cos  G  +  ya  sin  G 


loll385581* 

lo087U6588 

0.11377966 

-O0863I46I38 

0oll022l*8U 
Ool2277ii58 


1.07382217 

1„08887727 

Ooll607858 

=0o86l39582 

0.1111*0287 
Oo 11962687 


b^  =  tan  0     cos  G 

-0.19336597    -Oo21339059 
b£  =-tan  Q    sin  G 

-0.56513760      -Oo6098U827 
b->  =  -  cos  (CX  -00 

-0.9U751372  -Oo96501iil*5 

b^  =  sin  (OC'-Ct)  (R  sin  6  -Ba) 


=Oo3197l509 

0al)n?iMi30 

Oo 21092082 

=00069U7652 

0.02221269 


-0 «26219670 
Ool3637898 
0o21092082 

-0o07U51*l8it 
Oo0195UU62 


1.0181*8322 

1.09191551 

0.121091*1*1 

-0.8567271*9 

0.111*09361 
0.1162021*3 


-0„2563U782 


.0.6981*61*11* 


-0.990881*56 


-0.13U71371* 
Oo 12520007 
Oo 21092082 

-0.08572075 
0o  01151*776 


TABLE     19  (Continued)  62 

OBSERVATION  1       OBSERVATION  2  OBSERVATION  3 

bcj    =  sec^  £)(r  Cos  Q  -oca  cos  G  +  ya  sin  G) 

sec2  6    r  lo35677091  1.U7U50U6  1.&3S6636 

(R  cos  0  ~xa  cos  ®+ 

ya  sin  G)  0„11022U8U  0.H1U0287  0.111*09361 

b£  0.1U9&986  0.1^79080^  0.1772^199 


TABLE     20  63 

SOLUTION  OF  DIFFERENTIAL  EQUATIONS 
SECOND  APPROXIMATION 
OBSERVATION  1       OBSERVATION  2       OBSERVATION  3 


bin®- 

&6' 


Aa'  +  b*A6' 


-0.00009610 

-O.O0OO1739 

Oo 00000971 

-0.00000768 


0.00006317 

0.00001032 

-0.000051*29 

-0.00001*397 


-0.00000170 
-0„ 00000170 

0.0000^155 
0.00005138 


OBSERVATION  EQUATIONS 


dxa 

-1.055381*18 
-0.19336597 
-1.0336211*3 
-0.21339059 
-0.98561191* 
-0.25631*782 


al  dxa  +  *2  dy« 
bi  dxa  +  b2  dya 

d7a 

0.001*1*7722 
-0.56513760 

0.07381686 
-0.60981*827 

0.216915U* 
■.0.6981*61*11* 


=  a: 


=  a3  d«'  . 

+   b3  dza  =  b^  dCX  +  b£  d 


-0.91*751372 
-0.965011*1*5 
-0.990881*56 


Oo 00009610 
Oo 00000768 

-Oo 00006317 
0.00001*397 
0.00000170 

-0.00005138 


-1.05081086 
-I.706OO96I 
-0.95986771* 
-1.7882093U 
-0.768691*80 
-1.91*571*790 


NORMAL  EQUATIONS 


dx 


a 


dya 


dz 


a 


3.302280071 


0.12361*5699 
1.231668857 


0.61*3153009 
1.816085355 
2.810887350 


-Oo 000035501 
Oo 000000868 
0.000001203 
0.000000018 


l*.0690l*3279 

3.1711*00778 

5.270126916 

-0.0000331*12 


dxa 

dxa    3*302280071 

dya  -0.0371*1*2523 

dza  -0.191*760285 

0.000010750 

dx* 


REDUCED  EQUATIONS 
dya  dza 


0.12361*5699 

1.227039250 

4.1*601*29309 

=0.000001790 


0.61*3153009 
1.792001*081* 
Oo  0685311*01 
■Oo 000071631 

dz, 


-0.000035501 
0.000002197 
0„  000001*909 
0.000000017 


l*.0690l*3279 
3.01901*5531 
0.068536309 
0.000000016 


s.a  dya  —  a 

+  0.000020851  +  0.000102822     -0.000071631 


xa    0.3511*62753 
Axa -I- 0.000020851 

x0  + 0.3511*83501* 


ya        0.912606236  za      0.210920821* 

Aya  +  O0O00102822   ^za  -0.000071631 

y0  +0.912709058  z0  +  0.21081*9193 


6k 

V7  =  I 

where  n  is  the  velocity  of  the  mean  anomaly  ^  M.     The  resultant  "a"  from 
this  formula  has  the  value  1  ,,17190297  in  contrast  with  the  value 
1.17207018  which  was  computed  from  the  formula     a  =  q/l-e.     The  final 
result  of  geographic  coordinates  using  this  value  and  the  linear  method 
is  listed  in  the  summary. 

Second,  the  linear  method  was  used  again  in  which  the  eccentricity 
derived  from  the  mean  orbital  elements  was  used  instead  of  the  smooth 
value  for  e.     The  mean  eccentricity  for  the  time  of  observation  is 
O.10067li97o     The  resultant  coordinates  from  this  experiment  are  also 
listed  in  the  summary. 


Longitude 9 


Summary: 

Latitude ,  0 

Method 

Approximate 

12°  13°   23".21  N 

Differential  (1st) 

12°  19  03". 168 

Differential  (2nd) 

12°  lii"  UU"oll8 

Linear  (using  computed  k) 

12°  19  06».32l 

Linear  (with  mean  e) 

12°  19  U5"»276 

Actual 

12°  09   20".U 

68°  $9  07 ".86 

68°  $6*  LU"o336 

68°  26*  18".0U2 

68°  56"  11". 072 

68°  91 '  22". 278 

68°  20 «  1U».0 


60     THE  INTER-CONTINENTAL  TIE 
The  last  step  of  this  thesis  is  to  effect  the  *tiert  between  the 
two   stations,  the  geographic  coordinates  of  which  are  now  known «     This 
problem  can  be  considered  to  be  one  of  geometric  geodesy  or  what  is 
more  commonly  known  as  the  "determination  of  long  lines"  and  is  considered 
to  be  a  part  of  higher  geodesy<> 

Ho  Fo  Rainsford  states  there  are  two  main  problems  in  the  field  of 


long  lines 


6,  p.12 


a)  Given  the  geographical  coordinates  of  two  points, 
to  find  the  distance  and  azimuths  between  themj 
commonly  called  the  reverse  (or  inverse)  problem) 

b)  Given  the  geographical  coordinates  of  one  point 
and  the  distance  and  azimuth  from  it  to  another, 
to  find  the  coordinates  of  the  second  point) 
commonly  called  the  direct  problem,. 

Since  only  the  geographic  coordinates  of  the  two  points  are  known 

in  this  case,  the  reverse  problem  applies » 

The  following  data  are  given % 

Point  1  Point  2 

Curacao*  N.  W„  I0  Olifantsfontein,  South  Africa 

.,     0;   12°  1U'  W.118  N  0L  2$°  $V  3U".700  S 

X    68°  £6°  l8«.0l£  W  Xh2Q°  ^  £LM-1«>  E 

Problems     Find  the  forward  azimuth,    C^12>   at  Point  1;  find  the 
back  azimuth,   (X  21,  at  Point  2 5  and  determine  the  distance,  sp  between 
the  two  points o 

Since  the  reference  surface  chosen  for  this  paper  is  the  Inter- 
national Ellipsoid,  the  distance,  s,  between  the  points  will  be  one  of 
the  basic  geodetic  curves  called  the  geodesic     It  is  defined  as  [13,P«: 

a  line  on  any  mathematically  defined  geodetic  reference 
surface  that  defines  the  shortest  distance  between  any 
two  points  on  that  reference  surface 0     It  is  geometrically 
equivalent  to  a  great  circle  on  the  sphere     It  is  a  unique 
line  on  the  reference  ellipsoid  that  contains  the  normal 
to  the  ellipsoid  at  each  point <» 
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66 
Refer  to  Figure  ?•     It  can  be  seen  that  the  line,  s,  connecting  the 
two  points,  Pi  and  P2<,  contains  the  normal  to  the  ellipsoid  along  its 
entirety*     The  forward  and  back  azimuths,    (X  12  m^       21  respectively, 
are  shown  at  the  points <>     The  curvature  of  the  geodesic  in  the  diagram 
is  exaggerated  for  clarity  of  figure  only0 


Figure  7« 

It  should  be  noted  that  the  distance  difference  between  a  geodesic 
and  a  plane  curve  is  very  small  but  the  difference  in  the  azimuths  of 
the  two  would  be  appreciable .  As  the  distance  between  points  increases, 
the  difference  between  the  two  lines  (geodesic  and  plane  curve)  diverges 
rapidly* 

There  are  several  methods  available  in  geometric  geodesy  to  deter- 
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mine  the  geodesic  and  the  azimuths  but  because  this  particular  problem 
is  concerned  with  the  inverse  problem  and  the  distance  between  the  two 
points  is  approximately  6<,000  miles<>  then  a  choice  of  only  two  methods 
is  available.     These  are  the  Rainsford  Method  and  Sodano's  Fourth 
Method.     Another  method  that  solves  for  distances  greater  than  £,000 
miles  is  Rudoe°s  Method  but  since  this  procedure  treats  the  normal 
section  only,  it  is  not  applicable  for  this  problem.     The  method 
chosen  was  Sodano's  Fourth  which  is  a  non-iterative  solution.     The 
accuracy  is  of  the  order  of  iO" .1  in  azimuth  and  +  0<>5  meters  at 
6,000  miles o 


Ellipsoid 


Sphere 


Figure  8o 
Proceeding  with  the  solution 9  an  auxiliary  spherical  triangle  is 
employed.     In  our  problem,  the  westward  point,  Pi,  is  Curacao  and  the 
eastward  point,  P2  is  Olifantsfontein,  South  Africa.,     The  difference  in 
longitude,  L,  is  always  positive.  ^13,  p<>lp.     . 

The  first  step  is  to  find  the  parametric  latitudes,  pt  and  p29 
using  the  formula 

tan  B    -  (1  -  f)  tan  <£ 
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with  the  results 

fi  1    -     12°  07s   33"oS21 

fi  2  =-2$°  $V  Ol"ol89 
Using  the  parametric  latitudes,  the  approximate  spherical  distance,    (j  , 


is  found  from 


13sPoiil     s 


cos  6  =    sin/;,   sin  J3Z    +  cos  pi  cos /4  cos  L  =-0o 201730^06 


and 


sin  0/z     =   cos  fi,  cos  fii  sin2  ^    *  sin1    ,Q*     'z 


2_ 


=    0.77^1^8582 


The  resulting  value  for  6  is  101°  38 8  17"«336<, 

Now  A  must  be  found  using  the  formula    A  =  L  +  (   A  -  I»)  where)  13  $p.U2 
(  A  -L)  =  |fcl  +  Q  )  ^  *    (C+D)  sin  6   +     (F  +  E)  G 


„"  =  0°l8»20"o866 

■ 

where  P   -  20626Utt«80625o     The  value  for  Aiss 

L     =  97°  H?   O9"olli5 

(L~A)  =     0°  18"   20" .  866 

A    =   97°  29"  30" o Oil 

Then  the  forward  azimuth  0(i2  is  computed  from  the  formula 
[l3,  P0U2]   % 

cot  CX  12  =  [tan/3E    cos  fif  -  cos  A  sin/JJ/sinA  =  lHi°l6»06«o607 
The  back  azimuth  is     [13*  P«ii2  J  s 

cot «  21  [sin/?s  cos  A    -  coa  fit.  txifi,}  Sain  X    =  277°U98£6»o50U 
The  next  step  is  to  find  the  true  spherical  distance,    u  0,  using    A  = 
97°29J30"o011  instead  of  L    =  97°H909"oll45  in  the  same  formula  as 
before  s 

cos   0O  =  sin/?,    sinfi.L  +  cos  fii  cos fiz  cos  A  =  -0,2063 866912 
60  =  101°  5V  38".ii72 
The  final  step  is  to  find  the  geodesic  distance s  ss  from  the  formula 
[l3,   p.itf]    s 


69 

s   =b(A0   (90+    B0  sin    0o  cos  2  Cf  -  C0  sin  2  0  0  cos     1*   (/     ) 

s  =  11, 309,3U2„6228  meters     =    6, 102 „  14*21  nautical  miles <> 

Refer  to  the  following  tables  for  the  detailed  computations  of  the 

constants,  coefficients,  and  the  unknowns e 

Summarizing  the  results: 

Point  #1  Point  #2 

Curacao,  N.W.I.  Olifantsfontein,  South  Africa 

0,  12°  lli'W.118  N  Q>i&°  $1%  3U"o700  S 

\    68°  $6«l8"oOl£  W  X-u  28°  ^'   51" -100  E 


0(12    =  llit0  16 9  06»o607 
CX21   =  277°  U9«  56«o50U 


s    =  11,309S3U2.6228  meters, 
This  concludes  the  computations o 


TABLE     21  70 

LIST  OF  CONSTANTS 
FOR  USE  WITH  SCDANO'S  FOURTH  METHOD 

INTERNATIONAL  ELLIPSOID 


1/(X=    297.0 

a  =  6,378,388.0000  meters 

b  =  6,3^6,911»9U62  meters 

e  =  0.0819918898         e2  =  O0OO67226700 

e'  =  0.0822688896         e'  2  =  O0OO6768I702 


N  = 


/ 
e 


e*  +  e 


Kl=lfe4rr     =     237.2388918 

K2    =  16      — __  =  70^19.5HU5 

e2(l6e2N2  +  ef  * 

K*  =  I6e2  N2+e'  2      =  ^98652061*9 


2  e 


/  2 


^  =  iJl^e^  =  °^01081263° 

K*  =  l6el  ¥ _  =  0.799U593686 

l6e2N2  +  e/  2 

16  e2  *j2 — .  =  3.9865206U9 


K6  =:=1^7i 


TABIE  22  71 

COMPUTATION  OF  VARIOUS  TERMS 

A  =  cos'l     cos/^2.  sin  L 
sin  r\ 

cos  A     cos/Jl  0.8796086600 

sin  L  ^  0.99211*^600 

cos/}.    cosAsin  L      ,  O.8726998266 

(cos  yQ,  cos  yS^sin  iO/sin  g  0.8910182190 

B  =  A2   =   0.793913U666 

C    =    (cos  Q  -  B  cos   C7    wc* 
K3  U.9865206U9 

BcosJ  i  -0.1601^66007 

(cosg-  B  cos  8  )/  -0.01P.5739U99 

(cos  g  -B  cos  f)   MC3  =  -0.0083372662         =     C 


=    A 


=    D 


=     E 


D=  —  Ki;  sin  Bf    sin/3^ 

K,                      n  O.U010812630 

s4n#,     sin/3x  -0.091701^172 

-K^  sin/},    ainfl%  =            0.0367797603 

E  =    -K£  sinA    s±n.p% 

K*    a           r>  Oo799U593686 

sin  A  sin  Pi  =.0.091701^172 

-K£  sin  Qf  siny92  =             0.0733116370 

F  =  CK6 

C  -0.0083372662 

K6  3.986S206U9 

CK6  =          -0.0332366839          =     F 

G  =$2/sin(9 

a  101°38«17».336 

S  (radians)  1.7739203W* 

Q\         m  3.11*6793387 

Q/aind  =            3.2128U6105           =     G 


TABLE  23  72 

COMPUTATION 

PARAMETRIC  LATITUDE  &  APPROXIMATE  SPHERICAL  DISTANCE 

Parametric  Latitude  (  R) 

tan  R  =  (1-f )  tan  (/) 

f  Oo003367003U 

(1-f),  0.9966329966 

tan  Q)  0.21558157  -O.U8686087 

tan    Q  r,    0.21U8557061  -0.1*852216078 

/}     "         /J  =  12°  O7»33»o521  ^=-25°  53'01M„189 

APPROXIMATE  SPHERICAL  DISTANCE 

.    .  cos  0  —  sin/j,    sin  Zj^  +  cos  /?/ cos/\cos  L 

L  =AA  97°  ll«09»,li*5 

sin  L  0.99211*556 

cos  L  -0*12508862 

sin  /D  0.2100618662                                     -0. 1*3651*5281*9 

cos  P  0,9776880956                                     0.8996822851 

sin  n>  sin   O  -0.0917 015172 

cosfci  cos/^  0„8796086600 

cos/j/  cos  p   eos'L  =0.1100290331* 

cos  Q  -0.2017305506 


sin 


=     [cos^Q,    cos^  sin2  |  +  sin2     ^\  ^'  j 


2 

,o. 


L/2    T  ij8u3583U"o573 

sin    J  0.75002937 

sin2L?2  0.562514*0559 


Pirflk  -38°00«3litt.710 

^tr2-  ^ 

cos  ft,  cos/k  sin?  l/2     „-,  0.1*91*8186232 


{Rrk)/2  -l9°oo«  17"  .355 

sin  (  fo  -   o)/2  -0.32561*771 

sin2(/^-  ^0/2  0„1060i*61*310 


1&.  $/t 

R,   COSp^Si 

5  n  cos  />/  sin^  .£ 
fit       P^     n? 


cos  YD  cos  /•/  sin£  Jf  -f  sin 

A*-/? 


0.60086505U2 


sin     0/2  0.77515U8582 

@/2  50°1*9,08».668 

0     =  101°38817".336 


9  0.979UU1O57U 

cos  Q  -0.2017305506 


sin 


TABLE     2k 
COMPUTATION  OF     (  \  -  L),   (X  -^  and  OC  2l. 

(Kj  +    B   )    6   +     (C  +  D)  sin  (9+  (F+E)G 


(A-D 


Kl 
K2 
(Kl+8  ) 

»)  sin© 

:) 

(F  +  E)"     G 

r 


(C+D) 
(C+D) 

(F  +  E) 


K2 

237.2388918 
70£L9  .$!!}& 
238.03280^2666 
U22 .25123586 
O.O28l0;2U9Ul 
Oo0278U77U65 
0.01^007^31 
0.12875U6570 
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Ap« 


( A-  «  = 

A=L  +( 


U22.U078U83 

70519  05Ui;5 


(.8910182190)     (  P   )    =  0°  18 « 20". 866 
A-  L)  =  97°  ll«09n.lU5+0°l8«20M.866  =  97°  29»30".011 


cot  CX12   ~~    [tan/JjCOs  /i/-  cos  A  sinjj/j  /  sin  A 

cotCtip,    j  (~O.U8^22l6078)(.97768809^6)-(-.1303820U28)(0. 2100618662 )J 
~  0.^91^6383^1  ' 

cotOt12=   -O.Uii700709U5/0.991ii638331    =  -O.U508556737 
« 12  -     11U°  16'  06».607 

cotCX2i    =  fsin/?2  C0SA    -  cds J3i  tan  Bj\   /  sin  ^ 

eotCXo-.-  [(-*li36gU$28U9)(-0.1303820lt28)-(.89968228$l)(.21ii8g$706l)J 
<j.  0.9911*638331  " 

cotCX2i=      -0.13638U2066/0. 99Ha63 8331  -    -  0.137558U282 
CX2l=         277°  U9»  56".50U 


TABLE  25  Ik 

COMPUTATION  OF  TRUE  SPHERICAL  DISTANCE,  6 0. 

(USE  A  VICE  L) 

COMPUTATION  OF  GEODESIC 

cos  C0  =  sin  Pi   sin  /■>  2  -f  cos  Pi  cos  pz  cos  A 

=  (-.0917015172) +  (o8796086600)(-.1303820li28) 
=  (-.0917015172)  -(.11U68517UO) 
=  -o2063866912 

Q0      =  101°  51i«  38 

Q  0      =  1.778677025  radians 

cos  p0      a  I  cos  /j,    cos  pl  sin  A  J    /sin  (JQ 
cos  fl,  cos  fix  0o8796086600 

sin  A  0.9911*638331 

sin  g)0  0.978U705080 

cos  /3  0.8912891769  =  maximum  latitude 

P°  of  geodesic 

Ao     =  1  +  («'    2/U)  sin2y90     -  (  3A  )  8iaU^o 

e  2/h  OoOOl6920U26 

sin2/50  0.2056036030 

3e/ U/6U  0 .OO0O021U73 

sinU  /50  O.OU22728H16 

A0  1.0003li77993 

g0    =    (e  2A)  sin2^o    „(e'  U/l6)  sin^D 

e  2/t  0o00l6920U26 

sin2  po  0.2056036030 

e' Vl6  Oo 0000028630 

sink  fi0  OoOU22728hl6 

Ps9  0.0003U77691 

C0    =     •*       sinU/3 

e  Vl28  0.0000003579 

sinl*  0>o  O.OU22728U16 

C0  0.0000000151 

cos  2  O  _     2  sin  Pi  sinff^    -  cos  0 


sin* 

2  sinp,  OoU20123732U 

sin     pi  =O.U365It528ii9 
sin2  Po  0.2056036030 

cos    0O  0.2063866912 

cos     2ff  -0.6856357822 
2  0  033°17»08«.26l 

UG  266°3U516".522 

cos  h(f  -0.0598071587 
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TABLE  25  (Continued) 
s    =   b  (Aq   60  +  B0  sin    90  cos  2  (/ ~  C0  sin  2   0O  cos  1*  CT  ) 

A0  1.00031*77993 

0  o  1.778677025 

A0  0O  1.77929561*71 

B0  ^  Oo0003U77691 

s3ji   C70  0.978U7 05080 

cos  2(f  .  -0.6856357822 

B0sing     cos  2  0  -0.0002333091* 

C0  .  0„ 0000000151 

sin  2  jj0  -0.1*038865811* 

cos  k(f     r  ,  -0<,0598071587 

C0  sin  2    U0  cos  1*  U  OoOOOOOOOOOl* 

b  6,356,911.91*62 

s  11,309,31*2.6228  meters 

6,102.1*1*21  nautical  miles 

79027.5723  statute  miles 


7o  CONCLUSICN 

Upon  reviewing  the  results  of  the  determination  of  the  orbit  of 
Explorer  I,  it  can  be  stated  that9  in  general,  reliable  orbit  para- 
meters cannot  be  reduced  from  only  one  set  of  three  observations  of 
the  satellite o 

Based  on  the  results  of  the  calculations  of  i  and  u,  the  position 
of  the  orbital  plane  in  space  can  be  accurately  determined  from  a  few 
observations.  To  determine  the  size  and  shape  of  the  elliptical  orbit 
and  to  accurately  orient  the  orbit  in  its  plane  requires  a  series  of 
precisely  observed  and  reduced  coordinates*  As  Moulton  states  U*p.i4i  : 
the  difficulties  arise  in  finding  a5  e,  and  w0  As  mentioned  previous*- 
ly,  these  three  parameters  determine  the  size,  shape,  and  orientation 
of  the  orbit  in  its  plane.  To  determine  these  three  parameters  satis- 
factorily, the  whole  time  interval  of  the  three  observations  should 
be  divided  into  two  nearly  equal  parts  by  the  second  observation.  If 
this  is  so  then  0  ^   and  U3  will  be  nearly  equal.  With  the  time  inter- 
vals equal  or  nearly  so,  the  ratio  of  the  area  of  the  sector  to  the 
area  of  the  triangles  (see  Figure  5>)  will  be  more  reliable  and  the 
series  which  is  used  to  determine  the  auxiliary  parameter  p  will 
converge  more  rapidly.  Additional  terms  can  be  used  in  the  series 
to  determine  p. 

In  an  analysis  of  the  solution  for  ui,  ug,  and  u^  (derived  from 
p),  upon  which  (7^,  (f 2'   and  ^3  are  Dasedi>  the  reader  should  note  the 
value  of  each  u  determined  by  its  sine  and  cosines 

ul  u2  u3 

sine  value  21*8°  13 8 02". 933  2*3°  13»$0».317         2#°  3U82$».033 

cosine  value         21+8°  13 9 02". 7^6  2$3     13°50».30$        255°  li6'17"»8lfr 

difference  0«M77  0».012  11«52".781 
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There  is  a  close  agreement  between  the  value  of  the  angle  determined 
by  the  sine  and  that  determined  by  the  cosine  for  u^  and  Up  but  there 
is  an  11 « £2" » 781  disagreement  between  the  value  for  U3  derived  from 
its  sine  and  cosine »     The  value  used  for  U3  in  the  ensuing  calcula- 
tions was  an  average  of  these  two  which  is  2£?°  1*0»21".U2U.     This  could 
be  the  source  of  the  error  existing  between  the  computed  a,  ep  and  w 
and  those  calculated  from  the  82  observations «     I  believe  that  this 
can  be  attributed  to  the  uneven  time  interval  of  the  observations s 

TX        22h29m10Soll|0  ) 

T«         22  31  Oiio  U20  )       1  5l*so280 

T3    22  32  03  olOO  )   OPS8S.68O 

and  the  unrealistic  observed  declinations % 

Obso  #1        Obso  #2       Obso  #3 
hl°   22 »  S       51°  06 <  S      k9°   258  S 

It  seems  logical  to  assume  that  the  declination  should  either  increase 

or  decrease  progressively  during  the  short  intervals  of  the  observations, 

particularly  during  the  same  transit  of  the  satellite »  There  is  no 

such  progression  in  this  case0 

The  assumption  that  the  orbit  is  correct  is  another  source  of 
error »  The  accuracy  of  the  orbital  elements  predicted  for  the  time 
of  the  observations  from  the  unknown  station  depend  on  the  accuracy  of 
the  initial  values  of  the  parameters  determined  from  the  82  observations o 
The  standard  error  of  one  observation  of  +  £" -U3  undoubtedly  lessens 
the  reliability  of  these  values $   especially  the  mean  anomaly,  M,  which 
is  probably  the  crux  of  the  entire  problem,. 

The  assumption  that  the  observed  coordinates  of  the  satellite 
and  the  computed  orbital  elements  are  referred  to  the  same  mean  coordi- 


nate system  also  influences  the  result  <.  As  stated  in  95Po2 


on  the 
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orbital  elements 8 

the  reference  plane  is  defined  as  the  true  equator  of 
the  date0     The  origin  of  right  ascension  is  a  line 
shifted  from  the  mean  equinox  of  the  date  by  an  amount 
equal  to  the  precession  in  right  ascension  between 
19S0 pO  and  the  date» 

Given  below  is  a  formula  with  which  values  can  be  obtained  to  correct 

the  right  ascension  given  in  the  orbital  elements s  in  a  right  ascension 

referring  to  the  mean  equinox  of  the  date: 

Qp  (T)=  Q°  (DOl)  +  3°o$08  x  10=5  (Smithsonian  Day  -  33286) 
where  DOI  indicates  the  values  determined  by  the  Differential  Orbit 
Improvement  Program,,     This  information  was  not  available. 

It  is  apparent  that  the  methods  and  procedures  demonstrated 
herein  to  effect  the  tie  between  the  two  points  on  the  surface  of  the 
earth  can^  with  refinements 9  be  employed  in  practical  applications o 
The  success  of  this  application  depends  substantially  on  the  charac- 
teristics of  the  satellite  it self 0     In  facto  the  only  satellite  that 
would  be  practical  for  precise  results  would  be  a  geodetic  satellite  - 
one  that  would  have  the  following  characteristics s 

a»     Have  a  nearly  circular  orbit  or  one  that  has  a 

very  small  ellipticity  -  this  would  permit  reliable 
predictions  to  be  made  in  its  position  and  eventually 
an  ephemeris  could  be  made  available  3 

bo     Be  spherical  in  shape  -  so  that  its  motion  in  space 
would  not  affect  its  presentation  surface | 

Co     Have  a  perigee  of  approximately  1,000  miles  to 
minimize  the  effect  of  atmospheric  drag? 

do     Have  a  large  inclination  angle  so  that  it  would 
be  visible  from  most  latitudes! 

e<,     Have  a  flashing  beacon  that  would  be  visible  at 
its  long  range o 
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